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Abstract 

We provide a parametric representation for position-space Feynman integrals of a massive, 
self-interacting scalar field in deSitter spacetime, for an arbitrary graph. The expression is 
given as a multiple contour integral over a kernel whose structure is determined by the set of 
all trees (or forests) within the graph, and it belongs to a class of generalized hypergeometric 
functions. We argue from this representation that connected deSitter «-point vacuum corre- 
lation functions have exponential decay for large proper time-separation, and also decay for 
large spatial separation, to arbitrary orders in perturbation theory. Our results may be viewed 
as an analog of the so-called cosmic-no-hair theorem in the context of a quantized test scalar 
field. This work has significant overlap with a paper by Marolf and Morrison, which is being 
released simultaneously. 
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1 Introduction 

DeSitter space is a model for an expanding universe that is believed to describe both the earliest 
epoch of our universe, as well as the present epoch of accelerated expansion. It is therefore impor- 
tant to understand well the dynamics of classical as well as quantized fields on this spacetime. In 
quantum field theory, one wants to understand the behavior of correlation functions of the field(s) 
of the theory, such as the 1-point function (<|)a(*))*p> 2-point function (§A(x)§B(y))y, etc. where 
*F is some quantum state, and §a are some (possibly composite-) quantum fields. 

When the points x and y are a proper distance or proper time apart which is not larger than 
the Hubble radius, H , then one expects intuitively that the behavior of correlation functions 
in a generic state will not be substantially different from a "corresponding" state in Minkowski 
spacetime. This statement can be made somewhat more precise via the operator product expansion, 
which states that 

Mx)My)h ~ £cAfl(*,:y) (^cW)f ■ (i.i) 

c 

Each coefficient C^ B itself is known [|33l , at least in perturbation theory, to have an expansion in 
terms of the geodesic distance between x, y and curvature invariants at y. Thus, for distances or 
proper times much less than the Hubble radius (and masses greater than H), the coefficients will 
essentially be equal to those in Minkowski spacetime. As a consequence, for a state *P on deSitter 
such that the form factors ((|>c(.y))<p at point y are approximately equal to those of a corresponding 
state on Minkowski spacetime, the deSitter 2-point correlation function will not essentially differ 
from the Minkowski correlation function of that state near y. Thus, in this sense, one does not 
expect large departures of deSitter correlators from corresponding Minkowski correlators at small 
distances. 

For large distances, however, correlators could in principle show an unexpected behavior, and 
this intriguing possibility has been discussed by many authors in the literature, and for various 
theories, see e.g. B551 |27l HB HI E21 @71 HH, with partly contradicting claims. The simplest 
theory one can study is that of a free, minimally coupled Klein-Gordon field (|) with mass param- 
eter m 2 . When the mass parameter is strictly positive, then there exists a unique deSitter invari- 
ant Hadamard-state [Q3 H71 |9l, sometimes called "Bunch-Davies-" or also "Hartle-Hawking-" or 
"Euclidean-" state. For large proper time-like separations x, the 2-point function of this state de- 
cays as e~ M l T l for some M > depending on H and the mass, and a similar statement also applies 
to the connected higher «-point functions. Letting the Hilbert space vector representing the deSitter 
invariant state be |0), we may then also form a class of new states 

I^HE / fn{x h ...,X n )^X 1 )...^X n ) |0> (1.2) 

by applying a smeared product of operators to it, where the / n 's are smooth smearing functions 
vanishing outside some compact region. From the exponential decay of the n-point functions of 
|0), it then follows that e.g. 

(4>(*)>* = 0(e- M ^) (1.3) 

as the proper time separation x of x relative to any reference point goes to infinity. By the Reeh- 
Schlieder theorem ll54ll in curved spacetime, the class of such states l^) is dense in the Hilbert 
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space built upon |0). Therefore, this result states that the expectation value of (|) decays expo- 
nentially in time in generic states. A similar statement also holds for higher n -point correlation 
functions. An analogous result for the graviton field h^ v , with the appropriate value of M, could be 
viewed as a quantum analog of a "cosmic no-hair theorem". 

Because the graviton field is interacting with itself, it is of interest to ask whether a similar 
exponential decay property holds for the (connected) /7-point function in an interacting theory, e.g. 
when the self-interaction X§ 4 is turned on in massive Klein-Gordon theory. This question has been 
analyzed at 1-loop level in p4| by an explicit calculation, and it was found that the interaction does 
not change the exponential decay of the vacuum 2-point function. In the present paper, we consider 
this question for an arbitrary number of loops, for an arbitrary n-point function, and any dimension 
D > 2. We find that, again, the behavior does not change, so in this sense the quantum cosmic 
no hair theorem holds to all orders in renormalized perturbation theory in the analog system of a 
self-interacting test scalar field on deSitter. The decay constant M is the same in the interacting 
theory as in the free theory, i.e. it does not receive perturbative corrections. Our analysis is carried 
out at the level of individual Feynman graphs G, and the valence of the interaction vertices (= 4 
in ?u|) 4 -theory) is not important for our analysis. It therefore also applies to any other polynomial 
interaction. Of course, for a general polynomial interaction, and general dimension D, the theory is 
to be viewed then as an effective theory only, and the usual limitations concerning the interpretation 
apply. The same result has been obtained simultaneously, using a different method, by [[45 JU. 

In order to demonstrate exponential decay in time for the perturbative corrections to the corre- 
lation functions, we first develop a parametric representation for arbitrary renormalized Feynman 
integrals which is also of some interest in its own right. Our representation is based on the use of 
a Mellin-representation of the free propagator, which is combined with a version of the Schwinger 
parameter trick, and with some results from graph theory. The renormalization procedure is essen- 
tially that developed by l|29l l30l [TOl . specialized to deSitter spacetime and taking also advantage 
of many simplifications that occur due to the large symmetry group of this spacetime. Our fi- 
nal formula for the renormalized position space Feynman integrals Iq is eq. (14.691 ) as well as its 
Mellin-Barnes counterpart (15.961 ), which we repeat here: 



Tq is a kernel made from products of gamma-functions, and the quantities Z rs are deSitter space 
invariants related in a simple way to the signed squared geodesic distance between each pair of 
deSitter points X r and X s . There is one path of integration over a complex variable wp for each 
"forest" F in the underlying Feynman graph G, and the forests appearing in the exponent are 
precisely those connecting the external leg X r with X s . For example, for the following Feynman 
graph G in X§ 4 -theory with four external legs labeled Xi, . . . ,Xi, 

1 We are grateful to the authors for making available to us their manuscript before publication. 




(1.4) 
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-A', 



- A' 4 



an associated forest F in £4(2, 4) connecting X% with X4 can look like this: 

Xi 1 1 1 1 x 3 




As we shall explain, the forests are closely related to graph polynomials of the underlying Feyn- 
man graph ll56l . Such graph polynomials have also appeared previously in the study of Minkowski 
space Feynman integrals. However, our setup is different from that in Minkowski spacetime (see 
e.g. |H|42]|) due to the fact that (a) we work in position space, that (b) there are additional integrals 
of Mellin-Barnes type, and that (c) our graphs effectively involve an additional "virtual vertex", 
called "*" throughout this paper, as indicated in the above picture. The appearance of this vir- 
tual vertex is related to the fact that the loop integrations are over the deSitter hyperboloid, rather 
than flat space. To demonstrate the exponential decay of the connected correlation functions to 
arbitrary orders in perturbation theory, we use the above Mellin-Barnes representation and the ex- 
plicit knowledge of the poles of the meromorphic kernel Tq- Our Mellin-Barnes represenation also 
shows that Iq belongs to a class of special functions that have been called generalized //-functions, 
see appendix lAl 

In this work, we explicitly only treat massive scalar fields in the so-called "principal series." 
This assumption is made to simplify the discussion; massive scalar fields in the "complementary 
series" could presumably be treated by our method as well, and the only difference would be a 
correspondingly weaker exponential decay. However, our method would definitely fail for mass- 
less scalar fields, for which no deSitter invariant state exists even in the free theory (TJ. This case 
is a physically very important one, because interesting IR-effects are most often found in mass- 
less theories, and also because massless fields play a prominent role in various inflation scenarios. 
Some progress along these lines has been made e.g. by [1571 , but it appears unknown whether 
correlators in a generic, non-deSitter-invariant state decay exponentially in time or not. We also 
do not consider higher spin fields, such as the graviton field. For this field, a deSitter invariant 
propagator exists, but one has to deal with issues related to gauge invariance, see e.g. Il55l l26ll 
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for related discussions. We think that it would be very interesting to extend our analysis to such 
fields. Finally, it would be interesting to extend our results to non-perturbative constructions of the 
interacting Klein-Gordon quantum field theory on deSitter spacetime, which are available in two 
spacetime dimensions ll37ll . 

Notations and Conventions: A graph G is a collection of vertices VG = {1, . . . , V +E}, together 
with a set of lines EG = . . . , l r }, each of which connects a pair of distinct vertices. Our graphs 
have no "tadpoles". DeSitter points are typically denoted by capital letters such as X. d/u = \g\ 2 d D x 
is the invariant integration measure on D-dimensional deSitter or S D , depending on the context. 
*R(z) denote the imaginary- resp. real parts of a complex number. 

2 DeSitter space basics 

DeSitter spacetime can be characterized in a variety of different ways. Maybe the most natural 
description of this spacetime is as the embedded hyperboloid 

dS D = {X el D+l \x-x = -x^ + xf + ...+x^ = H- 2 } (2.5) 

in an ambient (D + 1) -dimensional Minkowski spacetime, with the induced metric. As is evident 
from this description, the isometry group of dSo is 0(D, 1), and it is hence a space of maximal 
symmetry H. DeSitter spacetime is a solution to the Einstein equations with a positive cosmo- 
logical constant A = j(D — l)(D — 2)H 2 , and with a vanishing Weyl tensor. It is entirely well- 
behaved from the point of view of causality — technically speaking, it is an example of a "globally 
hyperbolic spacetime". As a consequence, the initial value problem for an equation like the Klein- 
Gordon equation 

(V 2 -m 2 )(^ = (2.6) 

is well-posed, for any value of m 2 , including negative. Although deSitter space has the same 
number of Killing vector fields as Minkowski space, it has no analog of time translations, i.e. 
does not possess a Killing vector that is timelike everywhere. While this does not indicate any 
kind of causal pathology, it means that one cannot form a positive definite conserved Hamiltonian 
(energy) from the conserved stress tensor of the Klein-Gordon field, or similar other fields, 
although this is possible in certain sub-regions of deSitter spacetime, see below. In quantum field 
theory on deSitter spacetime, it also makes it impossible to define a reasonable global notion of 
particle with similar properties as in Minkowski spacetime. 

The deSitter hyperboloid admits various natural slicings and coordinate systems, some of 
which we consider in this paper. When we slice deSitter spacetime by a family of exponentially 
contracting, then expanding, round S 0-1 -spheres of constant Xq, the metric reads: 

ds 2 = -dT 2 +H~ 2 cosh 2 (Ht) d(D 2 D _ 1 . (2.7) 

2 This means that the number of Killing vectors, jD(D+ 1), is the maximum number that any (pseudo-) Riemannian 
manifold can have. 
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The coordinate x is the proper time of the geodesic curves of constant angle on the sphere. This 
coordinate system is useful e.g. to find the conformal diagram of deSitter spacetime. By making 
the coordinate transformation t = 2 arctan(e /f ' c ) G (—71, 7t) , the line element is brought into the form 

ds 2 = a 2 (-dt 2 + d® 2 D _ x ) = £T 2 ds 2 , (2.8) 

where ds 2 is the metric of the Einstein static universe, and where the conformal factor is Q.(t) = 
HcoshT l (H%). It is evident from this description that we can add a boundary to the deSitter 
manifold consisting of the two disjoint (D — l)-spheres labeled by J? ± = {t = ±7t}. While the 
deSitter metric is singular on this boundary, the "unphysical metric" ds 2 is smooth. For this reason, 
the boundaries are called "conformal boundaries" or "scri". 

It is sometimes of interest to consider subregions of deSitter spacetime as spacetimes in their 
own right. In cosmology, one is mostly interested in the "cosmological chart", which is the 
subregion of deSitter spacetime sliced by flat sections and which is covered by the coordinates 
(t,x) ER D defined by 

X = H 1 sinhHt - -He tH r 2 

2 

X\ = Q tH X\ 

Xd-i = e f xd-2 

X D = H^coshHt + ^He^r 2 . 

In this region, the metric takes the form 

ds 2 = -dt 2 + e 2Ht dx 2 (2.9) 

where dx 2 is the Euclidean flat metric on R D_1 . The cosmological chart covers the half {Xp — Xq < 
0} of dSj), and its boundary is sometimes called the (a) "cosmological horizon". The cosmological 
horizon is also equal to the boundary dJ + (i~) of the causal future of a point i of J 2 ^. The 
conformal diagram for the cosmological chart is: 



x = const. 




.D-l 



Another subregion of interest is the "static chart". It can be defined as the intersection of dSo 
with a wedge {X ±X\ > 0} in the ambient R D+1 . 
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static chart 




bifurcation 5 



■D-2 



The static chart is again a globally hyperbolic spacetime in its own right, and can also be defined 
as the intersection J + (i~) fl/~(z + ) of two points j -± £ <0 which are at the "same angle". It can 
be covered by the coordinate system (r|,r,x) defined for T| £ R, \r\ < H 1 ,x £ S D ~ 2 by 

X = (//" 2 -r 2 ) 1/2 sinh//ri 
Xi = (//" 2 -r 2 ) 1 / 2 cosh//ri 
X2 = rx\ 



X, 



D 



rx D - 2 



In those coordinates, the line element takes the form 



ds 1 



-f drf+r 1 dr 2 + r 2 d(£>l_ 2 , f=l-H 2 r 2 . 



(2.10) 



It can be seen from this form of the line element that, within this chart-but of course not in the 
full deSitter space-the metric is static, with timelike Killing field k = The corresponding flow 
T) 1 — y Y\-\- T defines a 1 -parameter group of isometries in the static chart, which correspond to a 
boost in the Xq-X[ plane in the ambient The boundary H = 9{+ U H- is formed from two 

intersecting cosmological horizons, and is an example of a "bifurcate Killing horizon", with Killing 
field k and with surface gravity K = H. This setup also occurs in the context of the black holes, or 
wedges in Minkowski spacetime, and is closely related with thermal effects in the corresponding 
quantum field theory on such backgrounds [|39ll22l[T5ll . 
Given two points X,X' £ dSp one can define 



H 2 X 



X' . 



(2.11) 



The quantity Z is sometimes called "point-pair invariant" because it is evidently invariant under 
the action of 0(D, 1). It is the analogue of the invariant distance squared in Minkowski spacetime. 
The causal relationships between points can be put in correspondence with values of Z, see the 
table [Hand the following conformal diagram: 
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Relationship between X and X' 



spacelike X' £ 7 ± (X) U 37 ± (X) 
lightlike X' G a7 ± (X) 
timelike X' G7 ± (Z) 



< 1 
= 1 
> 1 



Table 1: Shown here is the relation between Z and the causal relationship between X,X' '. 



Chart 


Metric 


Z 


Spherical 


-dx 2 + H- 2 cosh 2 (Hx) d& 2 D _ { 


cosh//(x-x / ) - (x -x'f cosh Hx' cosh Hx 


Cosmological 


-dt 2 + e 2Ht dx 2 


(e-^-e-^) 2 -^*-*') 2 


Static 


-fdx\ 2 + r l dr 2 + r 2 d(£> 2 D -2 


cosh[//(ri - if)] (f(r)f(r')) l 2 + H 2 rr' x ■ x' 


Hyperboloid 


-dX 2 + dX 2 + --- + dX 2 


H 2 X-X' 



Table 2: Shown here are the expressions for the point-pair invariant in various charts. 

J + (X) = future of X 



Z< -1 




north pole of 5 £ 



The relation to the signed squared geodesic distance a is given by 



cos(Hy/\a\) 
cosh(Zf-v/ja[) 



if points are spacelike, 
if points are timelike. 



(2.12) 



Table [2] records the expressions for the point-pair invariant Z in the various charts of deSitter 
spacetime described above. 

A spacetime metric whose corresponding unphysical metric = Cfig^v is smooth near J ;± 
is called "locally asymptotically deSitter'H. Such metrics have the property that they approach 
deSitter spacetime within a domain of fixed physical volume as that domain is moved towards 
J !± . The cosmic no hair theorem asserts that if we pick generic initial data for a metric on a 
spatial slice of the deSitter manifold and solve the Einstein equations with cosmological constant 
A, then such a metric will be locally asymptotically deSitter (away from any black holes that may 



3 A stronger notion of 
round sphere. 



asymptotically deSitter would be to require that the restrictions of ^ v to J !± is that of a 
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form). In this generality, the cosmic no hair theorem is still a conjecture, but a proof has been given 
when black holes are absent by [20] in D = 4, and by [|3l in higher dimensions. At the linearized 
level, the cosmic no-hair theorem says that, if h^ v is a perturbation of the deSitter metric, then we 
have (in a suitable gauge, see e.g. Il35l0 

h„ v = n% v , (2.13) 

where h^, is smooth (i.e. in particular finite) near Thus, in this sense, the cosmic no hair 
theorem is a statement about the growth/decay of perturbations near infinity. For other fields such 
as the Klein-Gordon field § on a deSitter background, a similar decay property holds. This is 
immediately evident for a conformally coupled field (m 2 = \D(D — 2)H 2 ); such a field has the 
behavior 

^ = a(^- 2 )/ 2 ^ hence ^ = 0(&- M ^) (2.14) 

where <}> is smooth near J*^, and M = H(D — 2)/2 for the conformally coupled field. For a general 
massive minimally coupled KG-field, M would be H(D — l)/2 — ^/H 2 (D- l) 2 /4-m 2 . 

3 QFT in general curved spacetimes and deSitter space 
3.1 General spacetimes 

The topic of this paper are quantum field theory correlators on deSitter spacetime in renormalized 
perturbation theory. However, to put things into a general perspective, we will begin by recalling 
some general key concepts regarding QFT in curved spacetimes. What we will say in this section 
is of a general nature and should apply to any (reasonable) QFT. For definiteness, the reader may 
think of the field theory described by the classical Lagrange density 

L = J/j[(V(^) 2 + fflV + ^ 4 ] (3.15) 

where d/j = \J—g d D x is the invariant integration element that is defined by the metric of a D- 
dimensional spacetime M. This is the theory that we will focus on in perturbation theory in the 
next sections, and our general remarks in this section are also going to apply to this theory, at least 
in the perturbative sense. 

A general curved spacetime does not have any isometries, so a quantum field theory on such a 
background will not have any of the features that are normally associated with these. In particular, 
one does not have conserved quantities like a hamiltonian, (angular) momenta etc., because their 
definition involves a time-like Killing field, translational (or rotational) Killing field, etc. For 
the same reason, one cannot define a vacuum- or ground state of the field, because its definition 
involves the hamiltonian. Without a vacuum (state of no particles), it is clearly also impossible 
to define in an unambiguous fashion particles, and in a spacetime without asymptotic regions 
approaching Minkowski spacetime, it is also clearly impossible to define an S-matrix, amplitudes, 
etc. 

While this means that one is forced to give up many familiar concepts when passing from 
Minkowski quantum field theory to a quantum field theory on a general curved spacetime, this 
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does not mean that quantum field theory on a curved spacetime as such is ill defined. In fact, in 
any spacetime with reasonable causal behavior (e.g. "globally hyperbolic"), we expect to be able 
to define a quantum field theory by (see 113311311 for details) 

1. The set of correlation functions (<|)(xi) • • •^(x n ))\p of the basic quantum field § (or more 
generally, any composite field) in a given state V I / . 

2. The operator product expansion [cf. eq. (11.11) 1 (OPE) between an arbitrary set of composite 
fields, whose coefficients are generally covariant. Because a quantum field theory will have 
many admissible states in general, the OPE can be used to determine that a given set of 
correlators belongs to the "same theory", but "different states", in that theory. 

An individual state should be thought of as being defined by the collection of its correlation func- 
tions. A quantum field theory will have infinitely many states, and their correlation functions will 
differ from each other substantially. But they will have the same OPE coefficients, as well as a 
number of other common properties which reflect general properties of any relativistic quantum 
field theory. The most important ones are that of "unitarity" or "positivity", and that of "causality". 
Positivity states that if A is any expression of the form 



with any finite number of smooth compactly supported testfunctions f n , then we should have 
(A*A)»p > 0, and we should also have the normalization (1)<j» = 1. The condition of causality 
is that, if x\ and xj are spacelike to each other, then any correlator containing [(|)(xi),(()(x2)] should 
vanish. The condition of causality is required in order that the theory has "propagation of infor- 
mation inside the light-cone'0. The condition of positivity is required for the probability interpre- 
tation of the state. This interpretation is as follows. Let A be a hermitian observable of the above 
form (13.161) . The probability distribution P^(a)da of this observable A in the state *F must satisfy, 
for all n, 



Under normal circumstances (i.e. when the moments /j„ = (A n )\p grow no faster than \/u n \ < 
est." n\), there indeed exists a unique probability distribution satisfying this formula for all n > 0, 
by the classical "Hamburger moment problem". The growth condition can be shown to be true to 
finite but arbitrary loop order in Xcp 4 -theory, and also non-perturbatively in flat spacetime in D = 2 
dimensions [43 J. The positivity requirement for the state implies that jU„ > for all n, and this in 
turn means that Pq>(a) > as required for a probability distribution. Furthermore, the normaliza- 
tion requirement of the state implies that the probability of measuring any value of a is equal to 1, 
again as required for a probability distribution. 

There are also, less obvious, regularity conditions that should be imposed on the correlation 
functions. These conditions state that, in a sense, the correlation functions are not just generic 

4 Note that this statement is entirely unconnected to the statement that the theory has "no tachyons". The latter 
statement is one about the stability of the theory, whereas the former is one about the propagation properties of 
signals. 



A= H $(*l)---4>(*n)/n(*l,-"> 



x n ) dn(x\)...dn{x n ) 



(3.16) 




(3.17) 
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distributions, but ones that are suitable "boundary values" of analytic functions^]. For a free Klein- 
Gordon field on deSitter, this condition is given below in sec. 13.31 see items 1), 2), and for the 
interacting X(|) 4 -theory in sec. 15.41 For Minkowski spacetime, such conditions follow from the 
"spectrum condition", i.e. the requirement that the Hamiltonian-or more properly, the energy mo- 
mentum operator-has non-negative spectrum. For a free field on Minkowski space, or a spacetime 
admitting a conserved non-negative Hamiltonian, another way of saying this is that the mode- 
functions defining the state are "positive frequency". While no conserved Hamiltonian exists in 
a generic curved spacetime, such regularity conditions can nevertheless be formulated in a very 
elegant and geometrical way using the mathematical machinery of microlocal analysis lfT2l . For 
the 2-point function, this condition is stated below in sec. 13.31 item 3). The microlocal condition 
may be viewed as a "tangent space" replacement of the usual spectrum condition and plays an 
important role in the renormalization theory on generic curved spacetimes 1121^1301 [roll - 
s' DeSitter spacetime 

If the spacetime admits any isometries, then particular states \P may be singled out by requiring 
that the correlators be invariant under the symmetries. Evidently, such conditions are most strin- 
gent in spacetimes with a maximally large isometry group, such as Minkowski or (anti-) deSitter 
spacetime. The invariance requirement is that (in the distributional sense) (§(gX{) ■ ■ ■ (j)(gX n )}\j( = 
(ty(Xi) ■ ■ -<])(X n ))<i< , for any element in the isometry group of the spacetime (or possibly the con- 
nected component of the identity if there are higher spin fields present), any X; and any n. In 
deSitter space, on which we focus from now orQ g G SO(D, 1)q. Because the deSitter group is 
non-compact, it is not clear that any such states will exist, and if so, if they will be unique. The 
free field theory example (X = 0) is instructive in this respect, and will be discussed in more detail 
below in sec. (|3.3I) . 

Because deSitter spacetime is an exponentially expanding spacetime, one might suspect that 
correlators of generic states in deSitter will decay exponentially for large time-like separation. 
More precisely, let us define the connected n-point functions as 

{Si,...,S r } jeSi jes r 

where the sum is over all partitions of the n-points into ordered disjoint subsets. As we show in 
sec. 15.41 in the context of perturbation theory, the correlation functions are analytic for configura- 
tions (X\,. . . ,X n ) such that no point is on the lightcone of another point. Consider a sequence of 
such configurations where one point X r -or more generally several points-go(es) to infinity in a 
time-like direction, i.e. in such a way that e " r ^- / Tjys — ^ 00 for nil s f, (hence |Tj-y| — \ Then 
we say that the state satisfies exponential clustering at large time-like distances if there exists a 
constant M > such that 

((()(Z 1 )...(t)(X n )}5 = 0(e- M ^l) ) (3.19) 

for any s ^ r. For a free Klein-Gordon field with m 2 > 0, exponential clustering indeed holds for 
the vacuum state and "most" other states, and the constant M is related to the mass and H, see 

5 If (M,g) is (real) analytic. 

6 In a theory like X(j) 4 which is additionally separately P, T-invariant, the invariance group is enlarged to 0(D, 1). 
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sec. 13.31 However, we note that it certainly not obvious that such a condition will still hold for 
general, i.e. interacting, QFT's. We will demonstrate this in sec. [6] in the context of perturbation 
theory for a massive Klein-Gordon field with self-interaction. 

As stated, the clustering property does not apply to configurations of points that are null-related, 
because for such configurations, the n-point functions are distributions. To obtain a sensible state- 
ment, one should smear the points out. One can easily formulate smeared versions of the ex- 
ponential clustering property. One such condition may e.g. be given as follows. Let A,B,C be 
observables as in eq. (13.161) . with testfunctions f n etc. supported in the static chart of deSitter 
spacetime, see eq. (12.101) . For T e R, let us denote by OCr(A) the time-shifted observable (13.161) . 
which is defined by replacing f n with /J, where we are setting /„ r ({r|,-, r,-,i,}) = f n ({T[i — 7\ r;, Jc;}). 
Then a natural alternative formulation of the time-like clustering property is that 

(Aa r (fl)C)£ = 0(e- M|r| ) as|r|^oo, (3.20) 

for any A, 5, C as in eq. (13.161) . 

This form of the time-like clustering property implies a quantum field theory analog of the 
cosmic no hair theorem, in the sense that all states approach the vacuum state at late times, at an 
exponential rate. Indeed, let (.)o be a deSitter invariant (i.e. vacuum-) state, which we assume 
satisfies the expontial clustering property as in eq. (13.20b . For some arbitrary but fixed A as in 
eq. (13.161) . let ( . )xp be the state whose n-point functions are given by^ 

WXO-.-^^ ^g^g^ . (3.21) 

Let Xj be in the static chart [cf. eq. (12.101) 1 with time coordinate T], in the static chart, and let Xj 
be the point which is time-translated by T into the future, i.e. it has time-coordinate T],- + T . Then 
it is easily seen from eq. (13.201) that, we have (in the sense of distributions in the static chart) 

Q{X?)--${x![)) v ^®{X 1 )-"${X n ))o asT^oo (3.22) 

at an exponential rate. In particular, for the 1-point function, we have (§(X T ))xp = 0(e~ M l r l), 
which is the analog of eq. (12.141) in the classical context. 

The clustering property can also be interpreted as saying that the probability distributions for 
two observables with large time-like separation become independent. Indeed, let A again be 
as in eq. (13.161) . with /„ supported within some bounded region in the static chart. Let ^ be a 
state with exponential clustering as in eq. (13.201) . and let A + = ar(A_) be the time-translated 
observable. Then from the clustering property, the n-th moment ((A + +A_) w )»p factorizes into the 
moments of A± up to an exponentially small correction. This implies that, for large times T, the 
probability distributions Pit for A± [cf. eq. (13.171 )1 become independent in the standard sense that 
the probability distribution of A + +A is given by the convolution of and . 

DeSitter spacetime is a real Lorentzian section of the complex deSitter spacetime dS^, defined 
by allowing complex vectors X in eq. (12.51) . 

dS% = {XeC D+l \X-X=H- 2 } . (3.23) 



7 In the GNS -representation of the state ( . )o, this state would be given by j^P) = A|0) up to normalization. 
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The complex deSitter manifold also contains the D-sphere S as a real Riemannian section. It is 
thus natural to ask what might be the relation between the QFT's associated with these two sections. 
In Minkowski spacetime, such a connection is well known and is described by the Osterwalder- 
Schrader (OS)-theorem [48J, which clarifies how to switch back and fourth between these QFT's. 
It essentially states that, under reasonable circumstances, the correlators are related by analytic 
continuation. There is also a corresponding theorem for the deSitter case [5J, and the main re- 
quirement for the correlators on the sphere is as in Minkowski spacetime a "reflection positivity" 
property. This reflection positivity property is the Euclidean counterpart of the Hilbert-space pos- 
itivity condition mentioned above. The analytic continuation from Euclidean deSitter space will 
also be important in this paper, but rather than appealing to the OS -theorem^, we will perform the 
analytic continuation "by hand" based on an explicit expression for the Euclidean correlators, see 
sec. [54] 

3.3 The free field on deSitter 

The free Klein-Gordon field provides a simple but instructive example to some of the general 
remarks made in the previous two subsections. It is defined classically by setting the interaction 
coupling X = in eq. (13.151) . The quantum field theory may be defined by its operator product 
expansion, but because it is such a simple theory, one can also directly give the relations in terms 
of commutators: 

1. The correlation functions of the free field in deSitter should satisfy 

(<|>(Xi) ■ ■ • [<b(Xi)^(X i+1 )} . ..<K*»)><p = iA(X h X i+1 ) (<|>(Xi) • ..tfXn)fo ( 3 - 24 ) 

for any state X V, where A is the commutator function on deSitter spacetime, equal to the 
advanced minus retarded fundamental solution for the Klein-Gordon operator (V 2 — m 2 ), 
see e.g. |fl9t 

2. We should have 

(V 2 - m\ (<|>(Xi) ■ • • 4>(X«) • ■ ■ tfX«))v = • (3.25) 

3. The state should satisfy positivity, {A*A)\p > for any A as in (|3.16l) . 

These conditions must be supplemented by a regularity condition in order to rule out states with 
pathological behavior such as the well-known "a-vacua'0 or "instantaneous ground states", for 
which e.g. the OPE does not hold. A common requirement which is essentially equivalent to the 
validity of an OPE in the case of a free field theory is that the state be "Hadamard". By this, one 
means that (i) the connected n-point functions are smooth (C°°) for n ^ 2, and that (ii) the 2-point 
function satisfies the first one of the following conditions: 

8 The OS-theorem works as it stands only for quantum field theories that are fully defined in the non-perturbative 
sense, whereas in this paper, we work in the perturbative setting. 

9 These deSitter invariant states 0] |47) have many pathological features such as infinite fluctuations for the stress 
tensor flTI . For a discussion of related problems, see e.g. l23l . The main deficiency of these states from the point of 
view of the present paper is that they are not boundary values of analytic functions, see 2.) below. 
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1. The 2-point function behaves like ((|)(Xi)(|)(X2))<p ~ £«„o^" + v log G e + smooth, where u n ,v 
are the Hadamard-deWitt coefficients, see e.g. Il46ll , and where a e = a + ite with ? the time- 
difference between the coordinates relative to an arbitrary foliation of deSitter by spacelike 
slices indicates the distributional boundary value (£ — > 0+). 

2. The 2-point function is the boundary value, in the distributional sense, of an analytic function 
defined on the neighborhood T 2 = {(Xi,X 2 ) G (dS%) 2 | 3(Xi -X 2 ) G V + }, where V ± is the 
future/past interior lightcone in IR D+1 , and where dSp = {X G C D+l | X - X = H 2 } is the 
complexified deSitter manifold. 

3. The analytic wave front set of the 2-point function is contained in {(Xi, /q;X2,&2) | 
k\j2 future/past directed} [[T21I321 . This version of the Hadamard property is a special case 
of the "microlocal spectrum condition" |[T2ll . 

The second and third conditions are essentially equivalent formulations of the first one [|52l[53l . 

It is also customary to restrict attention to so-called "Gaussian states", which are defined by 
demanding that the connected rc-point functions actually vanish for n^2. The n-point functions of 
such states are by definition determined by a version of "Wick's theorem" in terms of the 2-point 
function alone. It turns out that, for a massive Klein-Gordon field (m 2 > 0), there is a uniqu^^ 
Gaussian Hadamard state which is invariant under the full deSitter group 0(D, 1). This state, often 
called "Bunch-Davies state", has 2-point function Q3IUI 

<♦(*.)♦(*)).= ^m n ~ C) ll»J°~ l) * (-CD-l + c;D/2;l±^j, (3.26) 
where Z is the point-pair invariant, and where the dimensionless constant c is defined by 

D-\ l(D-l) 2 ^2 

+ V l -t l -hi <3 ' 27) 

As with any correlation function, this is distributional in nature, so it needs to be defined with 
some care as the boundary value (in the distributional sense) of an analytic function. The boundary 
prescription understood here is that explained in items 1.), 2.) above 1 ^ The Bunch-Davies state 
is sometimes referred to as a "vacuum state" because it has a comparable amount of symmetry 
than the ordinary free field vacuum state in flat Minkowski spacetime. However, one ought to 
be cautious when taking this analogy to literally. For example, in the static chart, where there 
is a natural notion of time-translation symmetry and a similar natural notion of Hamiltonian, the 
restriction of the Bunch-Davies state is actually a thermal state (a "KMS-state"), rather than a 
ground state, of this Hamiltonian 11391 . 

It turns out that the Bunch-Davies state can also be characterized as the unique analytic contin- 
uation (in the sense of boundary values of distributions) of the Euclidean Green's function defined 

10 If we drop the condition that the state be Hadamard, then there are in addition the a- vacua of (TJ 07), see also 
footnote|9] 

1 1 Namely we should replace X\ — > X\ + iee ,X2^X2, with e e V + , £ > 0, then smear out the distribution, and then 
let e ->• 0. 
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on the real Riemannian section S C dS^ to the real Lorentzian section dSp C dSp. This analytic 
continuation property is a manifestation of the OS -reconstruction theorem mentioned in the pre- 
vious section in this simple quantum field theory. Later, we will see that an analogous statement 
remains true for the n-point correlation functions of an interacting theory described by (13.151) . to 
all orders in X. An alternative "momentum space type-" expression for the Lorentzian 2-point 
function of the free KG-field which manifestly displays its positivity and analyticity properties 
has been provided by flU. Another alternative representation displaying explicitly the Hadamard 
property is given below in eq. (1A.120I) . yet a further representation is given below in eq. (14.311) . 

As in Minkowski space, there is an intimate relation to the representation theory of symmetry 
group of the spacetime, namely 0(D, 1) and the mass parameter, encoded here by c. We will not go 
into this topic here, but we note that the relation between the classification of the representations 
and c is (see e.g. 11581 ) 

1. c = +ip (pel), the principal series. 

2. c G [—(D — l)/2,0), the complementary series. 

3. c = 1, 2, . . . , the discrete series. 

The properties of the theory at large distances are different in each case; for positive m 2 , which is 
the case considered in this paper, only the first two cases can occur. For values m 2 < 0, no deSitter 
invariant vacuum state exists Q. 

The Bunch-Davies state satisfies the time-like clustering property (see above), which in the 
present free field theory reduces to the statement that for points X\ , X2 which are timelike separated 
by a large proper time x, 

((KXi)(KX2)>o = 0(e- M|x| ) (3.28) 

where M > 0. For the principal series M = H(D— 1 ) /2; this follows from a standard transformation 
formula for the hypergeometric function, or from the contour-integral representation given below 
in eq. (14.311) . For the complimentary series, the decay rate is only M = H(D — l)/2 — [H 2 (D — 
l) 2 /4-m 2 ] 1 /2 >0 . 

4 Interacting fields on dS-spacetime 
4.1 Basic setup 

In this section, we come to the main purpose of this paper, which is to give a prescription for 
defining the correlation functions (<))(Xi ) . . .^(Xn))^ for the interacting QFT corresponding to the 
classical Lagrangian (|3.15l) . We must face the following issues: 

1. We have to deal with any UV divergences in such a way that the renormalized correlators 
satisfy the desired general properties explained in sec. [3] 

2. We must specify which state ^ we would like to use. 

3. We must specify how to deal with any potential IR-divergence. 
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4. (Optional) We would like to check whether the time-like clustering holds. 

1) The renormalization problem has been dealt with for general curved spacetimes in |]29l[30ll . and 
we will apply this general procedure to the special case of deSitter spacetime in subsec. 14.21 2) 
For m > 0, it seems reasonable to use an 0(D, l)-invariant state, which in perturbation theory will 
be obtained from the Bunch-Davies state of the underlying free theory. 3) For this state, we can 
analytically continue the propagators from deSitter space dSp to the sphere S D . Working on the 
sphere will obviously eliminate any IR-divergences, and we will argue in sec. 15 .41 that the result 
can be continued back without introducing any new IR-divergences. 4) will be demonstrated in 
sec.[6l 

Following standard perturbation theory, our starting point will hence be the, as yet formal, 
expression 



OX 



AT 1 anal.cont. £ 1 U(Xi)^^) (x [ J 4 (Y) d^Y)) ) (4.29) 
xj-xtffl v> Vl \ V Js D J / Q 

The "expectation value" on the right side is intended to mean the Euclidean (analytically continued) 
Bunch-Davies state on the sphere in the free theory, and the expectation value on the left side that 
in the interacting theory on deSitter. N is an overall normalization factor making (l)o A, = 1> an d 
the fourth power (|) 4 in the free theory is defined via normal orderin g 12 l. Then we can (formally) 
apply Wick's theorem on the right side and get the standard expressior 
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Q(Xi)~-$(X E ))$fi = anal.cont. £ X v £ Sym(G) I G (X U . . . ,X E ) (4.30) 

Xi-+dS D y> Q G with V+E vertices 

for the connected correlators, where Iq is the contribution from a connected Feynman graph G with 
V interior vertices and E external legs, weighted by the appropriate symmetry factor (equal to the 
order of the automorphism group of the graph, | Aut(G) | _1 ). Formally, Iq is computed by writing 
down a factor of the Euclidean Green's function (ty(Xi)ty(Xj))o for each line (ij) E G, and then 
integrating the V "interaction vertices" over S D . These integrations lead to UV-divergences, and 
hence have to be defined in a more sophisticated way. In a general curved spacetime, this has been 
explained in 112911301 . and we now start explaining how the general procedure may be implemented 
concretely in the Euclidean deSitter spacetime, S D . For the A4 4 -theory, the graphs G have external 
vertices (Xi, . . . ,X E ) of valence 1, and interaction vertices (X E+ i,. . . ,Xe+v) of valence 4. But for 
most of what follows in the paper, the latter fact will not be so important for our discussion, so we 
may think of G as any (connected) graph with E external lines (without "tadpoles"). 



12 As explained in [29], normal ordering is actually not a renormalization prescription that is in general "local and 
covariant", nor "analytic in the coupling constants/metric". The correct definition of (|) 4 would be via the operator 
product expansion in the free field theory. This will differ from the normal ordering prescription by lower order terms, 
but can be compensated by making the re-definition (j) 4 — > (|) 4 + 6H 2 A (j) 2 + 3H 4 A 2 , where, with V|/ = T 1 /T and c as in 
the main text, A = (log// 2 //j 2 +\|/(— c) +V|/(c + 3))/16n 2 for some arbitrary /j of dimension [mass] (in D = 4). The 
advantage of normal ordering is that we can straightforwardly apply Wick's theorem. 

13 The normalization factor cancels as usual the contributions from bubble diagrams. 
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The first step is to bring the free Euclidean 2-point function (13.261) on the sphere S D of radius 
H , into a form that will be more suitable for our purposes. In expression (13.261 ), the two-point 
function is expressed in terms of (1 +Z)/2 = H 2 (X\ +X2) 2 /4, so the singularity of the Euclidean 
Green's function at X\ = X 2 is reflected by the singularity of the hypergeometric function at the 
value 1 of the argument. It is desirable to have a representation of this hypergeometric function 
which displays explicitly this singularity. A representation which is particularly useful for our 
purposes is given b^ A 



($(X 1 )$(X 2 ))o = Kc,d 



dz fH 2 (Xi-X 2 ) : 



c 2%i 



Here, 



K, 



c,D 



r{D-l+c + z)T(-c + z)T(-z) 
sin[7i(z+(D-2)/2)]r(D/2 + z) 

sm[n(c+(D-2)/2)}H D - 2 



(4.31) 



(4.32) 



(4it) D / 2 

the contour C is parallel to the imaginary axis near infinity, and it leaves the poles — (D — 2) /2 + No 
and No to the right, whereas it leaves the poles c — No and — (D — 1+ c) — No to the left. The contour 
C is visualized in the following figure, where we assume for simplicity that we have a principal 
series scalar field with c = — (D — l)/2 + z'p, and even D: 



C 



(D-l)/2 



-(P-2)fc 



5»(z) 



The above formula can be derived by noting that, for H 2 (X\ — X2) 2 /A < 1, the contour C can be 
deformed to a contour C' parallel to the real axis that encircles the (double when D is even) poles 
— (D — 2)/2 + No of the integrand. This integral can then be evaluated using the residue theorem, 
and it is seen to reproduce the expression (1A.120I) for the Green's function. We note that the line 
integral is easily absolutely convergent, as the integrand decays as e 2.tt;| | (h ere an( j i a t er we use 
the fact that |r(z)| ~ q-^)\/ 2 \3{z)\^- 1 ' 2 ). 

To obtain the Feynman integral Iq, we now need to multiply one factor of the Euclidean Green's 
function (14.311 ) for each line / of the graph G, and we need to integrate over the internal vertices 
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We obtained this formula for even D. The formula for all D (including odd) was obtained first by l45ll . 
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Xe+i, ■ ■ ■ ,Xe+v, while leaving the external points X\, . . . ,X E alone. We end up with an expression 
of the form 

I G (X h ...,X E ) = K^ D ( f] fdfKZ)) (U I ^)t 1 (X h ...,X E+v ) 

\i=E+l Js J \leEG JC Z7tz / 

x U2 2 ' n-c + zi)r(c + D-l+zi)r(-zi) 
L£ G sm[n(zi + (D-2)/2)]T(D/2 + zi) 

where L (= 2V +E /2 in (|) 4 -theory) is equal to the number of lines I in the graph G , and where t% 
is the product 

V+E 

t t {X h . . . ,X V+E ) = n [H 2 (X -Xj) 2 ] Zij . (4.34) 

In this expression, i, j is a pair of vertices, and for each such pair we have set 

Uj= £ zi, (4.35) 

l€G:l=(ij) 

i.e. Zij is the sum of all the parameters associated with lines that connect a given pair of vertices 
i, j = 1 , . . . , V+E. If there are no lines in the graph G connecting i, j, then the corresponding factor 
in eq. (14.341) is understood to be absent. 



4.2 Renormalization 

The ^/-integrations in the Feynman integral Iq [cf. eq. (14.331) 1 are not absolutely convergent for 
the values of zi that we need in the subsequent contour integrals for D > 2. The potential problems 
come from the divergence of the integrand t% near configurations of points such that X, = Xj for 
some i, j. The first task in this section will be to define these integrals over the Xi properly, by a 
method of extending distributions. This method is an application of the more general constructions 
in dTUl [301 |33l (based in turn on earlier work of lfT8ll ) to the special case of deSitter spacetime, 
taking into account many of its special features. 

To simplify the discussion, let us suppose that there are no external points, E = 0, so that Iq is 
simply a number. The basic type of quantity that we need to look at is (dropping H etc.) 

tc(f) := (fl fdfi{Xi)) f[ KX-XjfY'i f(X h ...,X v ) . (4.36) 

\i=l J S D J ij = i 

Here, we have put a superscript "G" to indicate the dependence on the graph G whose lines / 6 EG 
carry the labels zi G C. zij G C in (14.341) is given in terms of the parameters zi by eq. (|4.35l) . If 
there is no line in G between a pair i, j of vertices, then the corresponding factor is understood to 
be absent. / is some smooth test function on (S D ) . 

Of course, we will be interested in the case when / = 1 , but of course, this is a priori not well 
defined. But one notes that the integrand is actually perfectly smooth as long as Xi ^ Xj for all i, j, 
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so we can integrate it against any smooth weighting function f(Xi,...,Xy) whose support (i.e. the 
closure of all points where it is non-zero) is contained within the subset 



{(Xl,... ,X V ) G (S u ) v | X t ± Xj for all i ^ j} (4.37) 

of (S D ) v . Since we would like to take / = 1 , the problem at hand can be viewed as that of extending 
the distribution to, initially only defined on the above domain, to all of (S D ) V . The extension can 
then, by definition, be integrated against / = 1. In this view, UV-renormalization corresponds 
to extension of distributions, as originally proposed in [18J. Because we are not changing the 
integrand of to at the points (14.371 ) where it is already defined, it is intuitively clear — and can be 
rigorously proved [fT8l [TOl |29l — that this strategy for UV-renormalization corresponds to "local 
counterterms". 

To illustrate how the extension process works concretely, let us consider the following example 
in D = 4, and for V = 2 points: 

t{f) = / dn{Xddn{X 2 ) [{Xi -X 2 ) 2 f f{X u X 2 ) . (4.38) 

JS 4 xS 4 

As a distribution, this is defined a priori only for / which are nonzero only strictly within the 
set {X\ 7^ X 2 }, and the task is to extend it to a distribution called t R ("R" for "renormalized") on 
the entire product manifold (S 4 ) 2 , so that we can in particular set / = 1. In order to deal with 
the problematic configuration X\ = X 2 , let us introduce Riemannian normal coordinates around 
X 2 , i.e. we identify the tangent space Tx 2 S 4 with M 4 via a choice of orthonormal tetrad, and use 
the exponential map to identify a neighborhood of the origin in Tx 2 S 4 with a neighborhood of X 2 
in S 4 . The Riemann normal coordinates are denoted x = (x 1 , . . . ,x 4 ) and are defined in an open 
neighborhood of in IR 4 . As usual, the squared geodesic distance is given by a = x 2 . Then we 
have, using eq. (12.1 II) 

[(X.-X^Y = 2^[l - cos(HV?)Y = (x 2 Y (l - 2 £ 

oo 

= ^ aj H 2 j(x 2 Y +J , (4.39) 

7=0 

with easily computable coefficients aj depending analytically on z. Similarly, we can expand the 
integration measure as 

oo 

dniXi) = £ bj H 2j (x 2 ) j d 4 x , (4.40) 
with easily computable coefficients bj, and together this gives an expansion of the form 

oo 

t(X h X 2 ) dn{X{) = £ H 2j u z , 2j (x) d 4 x . (4.41) 

7=0 

The homogeneous distributions u Z)2 j of degree 2z + 2j are of the form 

»z,2j(f) = c } [d 4 x (x 2 Y +J f{x) , (4.42) 
Jm 4 
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for some easily computable numerical coefficients cj depending analytically on z. f is a testfunc- 
tion that a priori has to vanish in an open neighborhood of x = 0. But for sufficiently large values 
of the exponent, D\(2z + 2j) > —4, this distribution is automatically well-defined on all of IR 4 , in 
the sense that it has a unique extension to all of M 4 . So we need to worry only about the case when 
D\(2z + 2j) < —4, and in that case we simply proceed by analytically continuing the result from 
9\(2z + 2j)> —4. By lemma 6 and the following discussion of [|33l . we obtain in this way a family 
of extensions of u z ^ 2 j that is analytic in z G C, except at the simple poles 2z G —4 — No. Around 
those, u Zj 2j has the Laurent expansion 

N yi—2z—2 j—n 

Wto = I T^T d nJ (3 2 ) ( ^ 4)/2 8 4 W + regular, (4.43) 
^ 4 2z + 2j+n 

for some computable numerical constants dj >n , and an arbitrary constant M > with the dimension 
of mass. The regular part is determined by the condition that it is analytic in the half space D\(2z + 
2j)>-N. 

Because we have now extended the u z2 j for z such that 2z + 2j ^ — No, we also get a corre- 
sponding family of extensions of t to all of (S 4 ) 2 , which has simple poles for certain integer values 
of z. In fact, substituting the Laurent expansion for u Z)2 j into t gives, now in D dimensions: 

N {n-D)/2 M-2z-n H 2j 

t(X h X 2 ) = £ £ — — ©„ J , D (V 2 ,// 2 )6(X 1 ,X 2 )+regularpart 

n=D y=0 lZ + 11 

=: pp N [t(X h X 2 )}+ regular part. (4.44) 

Here, 'Dn.jp denotes a homogeneous polynomial in // 2 ,V 2 of degree (n — D)/2 — j with com- 
putable coefficients related to the d n j. The "regular part"by definition has the property that it has 
no poles for 9^(2z) > —N. The "pole part" pp N is the sum of delta functions on the D-sphere and 
its derivatives. In practice, we need to worry only about finitely many pole terms, i.e. finite N, so 
convergence is not an issue. The appearance of the 8-function on the sphere is clear because, if / 
vanishes near coincident points, then there is no divergence under the integral t(f), and hence no 
divergence at 2z = —D — No etc. either. Hence, a possible prescription for defining the extension 
of t(f) when D\(2z) > —N is to simply compute the integrals for sufficiently large 9t(z), then an- 
alytically continue, and then subtract the pole part (14.441) . after which one can simply set z to any 
desired value 9t(2z) > — N. In formulae, 

t R (X h X 2 ) := t(X h X 2 ) -pp N [t(X h X 2 )} . (4.45) 

Let us now continue this procedure to extend the distribution to of the form (14.341 ) for general 
V > 2. Again we are interested in the case / = 1. For this choice, the integral is singular, and 
we have to view it first as a distributional kernel that is, a priori, only defined on test-functions 
/ which are nonzero only strictly within the subset {X, ^ Xj for all i ^ j} of (S D ) V , where it is 
perfectly smooth. Again, we would like to extend this to the entire product manifold (S D ) V , so 
that we can integrate against / = 1 in particular, and so that the result is smooth (in fact analytic) 
in Zij in a suitable domain of C. We could again try to simply subtract the pole part as in the 
case of two vertices, but this would give the wrong answer, i.e. it would not be an extension but 
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a genuine modification of the distribution even at points where it is already defined. Furthermore, 
since the distribution t is now dependent on many complex variables zi, I G EG [related to the Zi/s 
in eq. (14.361 ) by eq. (14.351 )1, it is not so clear how the pole part etc. should be defined. 

The problem is that, in general, there are poles coming from any subgraph y C G (and sub- 
graphs thereof etc.), and these have to be disentangled. One way of proceeding is by induction. 
Let us associate a subgraph y(I) of G with any subset I c {1,. . . , V}, which is formed from the 
vertices in /, together with the set of all the lines in G joining them. If I\ , . . . , I r is a partition 
of {1, ... ,V} into pairwise disjoint subsets, then we let y(/i), . . . ,y(I r ) be the corresponding sub- 
graphs of G. Furthermore, we denote by G/ Ujy(Ij) the new graph with r vertices obtained by 
shrinking the vertices in each y(Ij) to a single vertex, and by G\ U/Y(//) the complement. For each 
I C { 1 , . . . , V }, let us furthermore choose a number Nj G Nq (to be further specified later except for 
Nj = if / has only one element), and 

zi= £ zi. (4.46) 

ZeEy(/) 

For y{I) the subgraph associated with / = {i, j}, the quantity zi coincides with zu defined above in 
eq. (14.351) . We define a finite subset Aq C C by induction in V, with the following properties: 

1. If G has only 2 vertices i, j, then Aq — {2zij}. 

2. Denoting P(V) the set of all partitions h,.. .,I r of {1,... ,V} excluding {1}, . . . , {V} and 
{1, . . . , V}, we have 

A G = {2 Z/ }U (J {A gMy( /,) -£>.}. (4.47) 

PeP(v),p=i>jij J 

For example, for the graph with V = 3 vertices and lines (12), (23), (13), we have Aq = {2(zi2 + 
Z23 + zn),2zi2 — N, 2z23 — N, 2zu — N}. where N = Ni is assumed to be the same for all / = 
{1,2},{1,3},{2,3}. 

Next, we define iteratively distributions t^, which should be thought of as "pre-renormalized" 
versions of tQ, where all poles from any subgraph have been subtracted. The distribution 
? G (Xi, . . . ,Xy) is characterized by the following properties: 

1. For G consisting only of V = 2 vertices {1,2}, we have ^(X^Xj) =tc(Xi,X2). 

2. If G has V vertices, is a well-defined distribution on (S D ) V \ {X\ = ■■■ = Xy}. 

3. has a scaling expansion of the form 

t G (X h . . . ,X V ) dfi(Xj) ~ £ £ M- 2z+ ?H 2k u 9 2 k (x u . . . ,x v -i) n d ° x J » ( 4 - 48 ) 

j=l peA c k=Q y=l 

where x\,. . . ,xy-i G (Tx v S D ) v ~ l are the Riemannian normal coordinates, identified with 
vectors in a neighborhood of the origin in MP, of X\, . . . ,Xy~i for fixed Xy. Each Un,2j is a 
homogeneous distribution of degree p + 2j. 
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4. If / has its support in (S D ) V \ {X\ = ■ ■■ = Xy}, then t G (f) is analytic in the variables zi, I G 
EG subject to the condition that, for any proper subset / C {1, . . . , V} we have 

A y(/) C {w G C | «K(w) > -iV/} . (4.49) 

5. ?g agrees with tc for configurations (X\ , . . . such that Xi ^ for all i ^ j. 

From items 2),3), one can show by the same arguments as in sec. 3.3 of 11331 that t G has a unique 
extension to all of (S D ) V . This extension is analytic in the zi G C, I G EG subject to condition (14.491 ) 
in 4) for all I, except possibly for those configurations zi G C such that Aq PI (—No) 7^ 0, where it can 
have poles at p = p(z) G Ag H (—No). The residues at these poles are proportional to derivatives of 
the delta-function 8(Xi , . . . ,Xy). These correspond to the "superficial divergence(s)" of the graph 
G, and we have in fact 

N n-{V-l)D M -2 z -n H 2j 

t G (X h ...,X v ) = £ £ £ — - !P nJijD| p8(Xi, . . . ,X V ) +regular part 

peA c «=(V-1)D ;=0 lPW+ n J J " H 

=: p PiV [?g(Zi,...,Zv)]+ regular part. (4.50) 

Here l D n ,j,D,p is a 0(D + l)-invariant partial differential operator of degree n — 2j — D(V — 1) on 
(S D ) V which is homogeneous in V,//, which may depend analytically on the z\ G C. The regular 
part has, by definition, no poles for any configuration of zi G C such that £H(Aym) > — Ni for any 
subset/, including now / = {1, ... ,V} itself (i.e. y(7) = G), where Nu y} '-=N. The extension 
to is then defined as this regular part, i.e. by performing the "MS-subtraction" 

t G -(X h ...,X v ):=t G (X h ...,X v )-pp N [t G (X u ...,X v )} . (4.51) 

What is left is to give a recursive definition for the pre-extension. The recursion starts with V = 2, 
for which the pre-extension is known, and then increases V. It is in fact uniquely determined in 
terms of l)-5) above and is given aa 15 l: 

t G (X h . . . ,X V ) := t G (X h . . . ,X V ) - (4.52) 

PeF{v),p=Ujij j 

This formula requires some comment. First, since the partition I = {1,...,V} is by definition 
excluded in the sum over partitions, the unknown term pp W/ [t G (Xk) ^ e /] does not appear on the right 
side. Furthermore, the partition P : {1}, . . . , {V} is excluded as well, so at least one Ij in each 
partition P ..,I r has more than one element. Therefore, the unknown term t G (Xi , . . . ,Xy) also 
does not appear in the right side. Hence this formula defines t G in terms of quantities that are 
already known inductively. 

15 When we iterate the extension procedure that we have sketched to general V, one ends up with a compact formula 
which is called "Zimmermann forest formula", because a similar formula was derived for momentum space Feynman 
integrals in Minkowski space in a famous paper by Zimmermann (59]. In position space, a formula of this kind has 

been obtained by HUE] US. 



24 



Because PPjVy.t^/.^^O&ei)] is proportional to derivatives of the 8-function with \Ij\ entries, this 
effectively means that the vertices from Ij in t^ u ^^(Xi,. . . ,Xy) are "merged". In fact, if there 
were no derivatives at all, we could replace this by tQ, u ,j^Xj l ,...,Xj r ), where Gj Ujj(Ij) is 

the graph with r < V vertices obtained by shrinking all graphs y(Ij) to a single point. If there are 
derivatives, the story is similar, except that the derivatives will first act on the points X^, k E Ij for 
each j = 1, . . . , r. This will produce again a distribution of the form , J JXj 1 ,. . . ,Xj r ) which 
has slightly shifted zi due to the action of the invariant operators <D j. One may be shown that the 
iterative definition satisfies l)-5) by induction — the proof of this is similar to that given in 11301 . 

As we have described above, the analyticity domain of as a function of the parameters 
Zi, I E EG, is given by those zi for which 

£H(Ay(/)) > —N{ for all /C {1,...,V} . (4.53) 

Thus, by taking all the Ni sufficiently large, i.e. by subtracting enough of the poles in the MS- 
subtraction step at each order, we can achieve that has an analyticity domain in the zi EC that 
is as large as we desire. 

In summary, we have outlined the proof of the following: 

Theorem 1. The distribution tc(X\, ...,Xy), initially only defined as a distribution non-coinciding 
configurations of points in S D , i.e. X ( - ^ Xj for all i ^ j has an extension to a distribution defined 
for all configurations (Xi, ...,Xy) E (S D ) V . For any arbitrary but fixed W G 1, we can define this 
extension so that it is analytic for zi E {w E C | *R(w) > —N} for all / E EG. 

Remark: 1) Although this result has been stated for graphs G without external legs (E — 0), it is 
easy to generalize the above argument to the case where external legs X\,...,Xe are present, as 
long as these do not coincide, Xi ^ Xj for all < i ^ j <E. 

2) In our formula for Iq, we need to integrate the complex parameters zi associated with the lines 
/ of the Feynman graph G over certain contours, and the location of these contours determines our 
choice of the numbers Nj. The contours are located for £H(z/) between — (D — 1)/2 and — (D — 2)/2 
for a principal series scalar field. Our extension to must be analytic for such zi- and we must make 
a choice of the Ni such that this domain contains the contours of interest. 

Given that we have now an extension defined in terms of to via subsequent MS-subtractions, 
we may ask whether we can go back and write to in terms of t^. For this, we just have to undo 
the subtractions at each step. This is a problem of basically combinatorical type, the solution to 
which is known in the literature on renormalization theory as "Zimmermann forest formula", for a 
derivation of this formula in position space, see [41 J. To state this formula in reasonably compact 
form, let P : I\ U • • • UI r = {1, . . . , V} be a partition of the set of integration vertices. Then define 
the inverse of the "MS-subtraction operator" by 

-7> ( g( Xl ,...,X v ) :=^ I) (X 1 ,...,X V ) pip Ko^Wl «J« (4J4) 
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A "forest" 4> over { 1 , . . . , V } is, in this context, a nested system of partitions, i.e. the first element 
P of the forest is just a partition as above, the second layer is a partition Pj further subdividing each 
set /, in the partition P, etc. The formula expressing t G in terms of is: 

t G {X h ...,X v )= £ (H{-T P )\t%{X l ,...,Xv). (4.55) 

forests <& V/>e<5 / 

In the product of the subtraction operators, the factors are ordered in such a way that the coarsest 
partition stands to the left. The forests are over partitions of the V interaction vertices of the graph 
G into subsets, which together with their lines can be thought of as sub-graphs of G. 

The right side of this is well-defined as a distribution on all of (S D ) V apart from the poles 
z/'s for which some p = p(z) £ Ay(n becomes a negative natural number for some / C {1, . . . , V}. 
Using the construction of the sets Ay, it is seen that this situation can occur at most when the z/'s 
are "in resonance" 

£ n izi £ Z for some £ Z. (4.56) 

leEG 

So we conclude: 

Theorem 2. The distribution tc(Xi, ...,Xy), initially only defined as a distribution non-coinciding 
configurations of points in S D , i.e. X, ^ Xj for all i ^ j has an extension to a distribution to all 
configurations, for z/'s such that the absence of singularities condition (14.561) is fulfilled. This 
extension is an analytic continuation in the z/'s of the distribution to in eq. (14.361) defined automat- 
ically for sufficiently large £H(z/)'s. 

Thus, to properly define Iq in formula (|4.33l) . we have two possibilities: We could use the 
extensions t^ provided by thm. \T\ Alternatively, if the z/-integrations are chosen parallel to the 
imaginary axis, so that the absence of resonance condition 

£ n$K{zi) £ Z for any n { £ Z, (4.57) 

ZeEG 

holds [in addition to the earlier restriction on the contour C described in the picture after eq. (14.311) 1 , 
then we may also use the extension tc supplied by thm. [21 The difference between these two 
prescriptions for defining Iq differs from the contribution of the pole terms. The pole terms contain 
8-functions and their derivatives in the X,-'s, see e.g. eq. (14.501) . and hence would give rise to 
additional finite residues in the z/-integrations proportional to 8-functions and their derivatives. 
Hence, by the "main theorem of perturbative renormalization" 031 [T3Tl . the results for Iq would 
differ from each other by terms that can be absorbed as local counterterms in the Lagrangian (I3.15I ). 

The upshot of the entire discussion is hence that we may simply use the formula (14.331) for 
the renormalized Iq, provided that the d/j(Xj) integrations are performed first, and the dzi inte- 
grations afterwards, where the analytic continuation of to in the z/'s is understood, and where it 
is understood that the contours of the z; -integrations satisfy the absence of resonance condition. 
That condition is satisfied if 1 and Zrl(zi),l £ EG are linearly independent over Q. This automat- 
ically incorporates a renormalization prescription compatible with adding/subtracting finite local 
countertermj^l . 

16 In D = 4 dimensions, these local counterterms would be of the standard form A(V(|)) 2 +B H 2 § 2 +C (j) 4 . Each 
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4.3 Parametric representation 

We next derive a parametric representation of Iq, cf. eq. (14.331) . As explained in the previous 
subsection, our formula Iq has to be understood in the sense that we first perform the integral 

,X V+E ) (4.58) 

for appropriate values ofz such that it is well-defined, then analytically continue, and then integrate 
over z along certain contours parallel to the imaginary axis subject to the absence of resonances 
condition (|4.56l) . In this section, we explain how one can do this in practice using the familiar 
trick of Schwinger parameters. The derivation is only formal, and the resulting integral over the 
parameters is not well-defined a priori. But it can be given sense via analytic continuation in z and 
a technique called "sector decomposition" [@] |24l [25J |7]. Another, this time rigorously derived, 
parametric representation based on the technique of Mellin-Barnes integrals will be provided below 
in subsec. |5~2| 

The first step is to replace the integrations over a D-sphere rather to ones Euclidean space. This 
can be achieved using the identity, 

dfi(X) = y 8(H 2 X 2 - 1) d D+l X , (4.59) 

together with 

b(H 2 X 2 - 1) = 6(log// 2 X 2 ) = f—(H 2 X 2 Y (4.60) 

J 2%i 

where the integration path is running parallel to the imaginary z-axis, so that the second expression 
is up to a prefactor just the ordinary Fourier transform of the delta-function. We would like to 
substitute this identity for each integration measure d/j(Xi) , i = E + 1 , . . . , V + E in the master 
integral, with new complex integration variables which we call z*;. We get, up to a prefactor 

(H/2) v , 

' V+E 



| ] / dn(Xi) t 2 (X h ...,X v+E ) (4.61) 

We next introduce Schwinger parameters for each factor under the integral (14.611) [cf. eq. (14.341) 1 
and use the standard formula 

[H 2 (X, - Xj) 2 } z = ——- ^ da a- l - z exp[-a// 2 (X, -Xj) 2 } . (4.62) 



A,B, C is a formal power series in X, whose coefficients can depend on c and \ogH 2 /M 2 . It is important to note that our 
prescription described here is not one that would in general have an "analytic dependence on the metric", in the sense 
of ||32]| . However, by adding appropriate finite local counterterms, we can cancel out this non-analytic dependence 
(in the present context, it would arise through powers of \ogH 2 /M 2 .) This follows because we know from [ 33 1 that 
a prescription with an analytic dependence of the metric exists, and that it will differ from the given one by local 
counterterms. 
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Inserting this into in all places clearly reduces the integrand to a Gaussian. We have one 
Schwinger parameter for each complex variable Zij or z*i, which we denote accordingly by ay 
or a,*. The collection of all Schwinger parameters is denoted a. The exponential in the Gaussian 
is minus 

V+E V+E 

£ a ij H 2 (X i -X j ) 2 + £ (Xj*H 2 xj (4.63) 

!J=1 j=E+l 
V+E V+E 

=■ £ QijtfXi-Xj + l £ Bj HXj + C. 

U=E+l j=E+l 

The matrix, vector and scalar quantities Q, 5, C are defined by the last equation. They are functions 
of the parameters a and B is additionally a function of the points X = (X\, . . . that are not 
integrated in Iq, i.e. the external legs. Explicitly, we have 

{Y.k^jk if i = j, sum over k connected to j, 
-OLij if i ^ j and i connected to j, i, j G {E + 1 , . . . , V +E} , (4.64) 

otherwise, 

as well as 

E 

Bi= £ ctijHXj, C = £ £ ay. (4.65) 

je{l,...,E}:(ij)eG i=lje{l,...,E}:(ij)eG 

These integration parameters can be decomposed conveniently as R ■ a, where R is a non-negative 
radial coordinate, and where a is now allowed to range through the standard simplex A, defined by 
the conditions a ;j > and Y,i = 1 where the sum is over those i,j such that either i = *,j = 
E + 1 , . . . ,E + V, or over i, j such that there is a line in G connecting i, j. 

Carrying out the d D+l Xi and the dR integrations is now straightforward, with the final result 
up to a constant not depending on z: 

s r(-z-(D + l)V/2) 
t ? (X h ...,X v+E ) = 1 ^ '-^ (4.66) 




Yin 



-Zij) 



1-J 



I X I 1 \ /v . . I 



Here and in the following, the quantity z is defined as 

z = £zy (4.67) 

«>;' 

and here, as well as in all the sums/products in eq. (14.691) . a sum £ ; ; (or a product) runs over i, j 
such that either / = *, j = E + 1 , . . . , E + V, or over /, j such that there is a line in G connecting i, j. 
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In the notation from the previous subsection, z = zi,I = {1, • • • , V}. The polynomials 11 , J depend 
on the graph G and are defined byS 

11 = det Q , 7 = det Q (C-^iQ'^ij Bi 'Bj) ■ (4-68) 

We will give another prescription how to obtain the polynomials 11 , y in the next section, where 
we also discuss many of the special properties that these polynomials have. Here we only note 
that the polynomials 11 , J are homogeneous in the variables oc^, and that each monomial of 11 has 
coefficient +1, whereas each monomial in J has coefficient +1 or coefficient (1 — Zyy), where Z^/ 
are the point-pair invariants formed from the external points X\,.. .,Xe- Hence, in the Euclidean 
domain where Z# G [—1,1), the coefficients of both 11, J are positive, and singularities in the 
integration over the parameters a in eq. (14.691 ) can therefore only arise from the boundary of the 
integration region A when at least of the a,j = 0. Because of these singularities, the integral 
eq. (14.661) is actually ill defined, and this is a consequence of the fact that our manipulations given 
above were in part only formal. 

We ignore this important point however for the moment, and we insert our expression (14.661) 
for the master integral into Iq, leading to 



^,'})^oi(ni'|f)<.-i« y ) 



X 



r(- z -(D+i)v/2) n 2 _ 2z , r((D - l)/2 + ip +zi)r((D - l)/2 - ip+zdn- Zl ) 
n r (-^) Mo sm[n(zi + (D-2)/2)]r(D/2 + zi) 



hi 

7 (D+l)V/2+z 
u {D+\){V+\)/2+z 



n a i/ ij ■ ( 4 - 69 ) 



This is the desired parametric representation for Iq. z stands for zi, I G EG, orz*i,« G {E + l,. . . ,E + 
V}, and /- means a multiple contour integral. For I E G the contours are parallel to the imaginary 
axis with 

-(D-l)/2<<H(z/)<-(D-l)/2 + e (4.70) 

for an arbitrarily small £ > . The remaining contours can be chosen at this stage to be any contour 
parallel to the imaginary axis. Kg is the constant 



Ka=H L(D-2)-VD 



coshTCp" 11 

.(471)^/2 



2 -V n V(D+l)/2 ^ (4J1) 



with L = | EG | the number of lines in the graph G. For the graphs in ?u|) 4 -theory, we have L = 
2V + E/2, so the dimensionful term is H V ( D ~ 4 " I+E ^~ , which in D = 4 reduces to H E . 

It is possible to make sense of the singular integrals over the Schwinger parameters via a 
method which has been called "iterated sector-decomposition" in the literature. It is normally 



17 Polynomials similar to 11 , J also appear in the context of loop integrals on Minkowski space, and are sometimes 
called "Symanzik polynomials" there, see e.g. Il36ll . However, we note that our polynomials are not identical to these, 
because we work on deSitter space (actually S D in this section), and because we are in position space. 
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carried out in the context of momentum space loop-integrals in flat space QFT. However, ex- 
pression (14.661) is of a form to which this formalism applies. Later in subsec. [5T2l we will present 
another parametric representation which avoids these problems, so will not discuss these ingenious 
constructions here and refer the reader to the literature 1141 12411251 171. 

We close this section by pointing out another more geometrical (but still formal-) way of writ- 
ing Iq. The a-integration can be viewed as an integration on the simplex A, where the integration 
measure is simply 



[da] 



1 



a, 



ij 



f\da 



i j ■ 



(4.72) 



Alternatively, let us give the simplex A an orientation by picking an order e\ < • • • < e r of the edges 
in EG*, and let Q.? be the differential form in a defined by 



j (D+\)V/2+z 
u {D+\){V+\)/2+z 



n 

i=l 



a? 



-I) a ek da ei A 



Ada ek ■ ■ ■ Ada e 



(4.73) 



Then a straightforward calculation based on the homogeneity of 11 , ¥ shows that 



d£l? = . 



(4.74) 



Let A be any smooth manifold such that dA = dA, and such that A is homotopic to A. Then the 
parametric form of the Feynman integral can be written as 



I G ({Z lj })=K G 

>< ri2 

leEG 



T(-z-(D+l)V/2) 



z J A 



II,,n :,,) 

2zi r((D-l)/2 + ip+zi)r((D-l)/2-ip + zi)T(-zi) 
sin[7i(z z + (D — 2)/2)]r(D/2 + zi) 



(4.75) 



5 Analytic continuation and Mellin-Barnes represenation 

5.1 The graph polynomials 11 , 1 

In this subsection, we describe begin describing the properties of the graph polynomials^ Uq-, 7g 
introduced in the previous section, which appear in our parametric formula (14.691) . Some of the 
properties of these polynomials will be used in the next section. 

Let G be a Feynman graph, with internal vertices i = E + 1 , . . . , V +E, and E external lines. In 
the following, we consider these external lines as part of the graph, and we let the external vertices 
on which they end be labeled by i = 1,. .. ,E, and each such vertex is associated with a X, 6 S D . 
A pair i, j of internal vertices can be connected by none or several lines. If they are connected by 
at least one line, then we associate a Schwinger parameter a i; 6 R to it. Since all that matters in 
the following is the association of the graph G with the Schwinger parameters, it is convenient to 
replace G with a graph G* such that 



In this section, we put a subscript "G" in order to indicate the dependence on the graph. 
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• G* has the same vertices as G, and an additional vertex called 

• Two vertices ij'G {1, . . . , V + E} in G* are connected by a single line if these vertices are 
connected by at least one line in G, and they are not connected otherwise. 

• Each of the vertices i = E + 1 , . . . ,E + V is connected to * with a line. 

With each of the last type of line, we associate a Schwinger parameter a,*. An example of a graph 
G and the corresponding graph G* is given in the following pictures: 





Following standard constructions in graph theory 11561 . let us define the Laplacian of the graph 
G* to be the V + E + 1 -dimensional square matrix (for i, j = * , 1 , . . . , V + E) 

OLjk if i = j, sum over k connected to j, 
—OLij if i 7^ j and i connected to j, (5.76) 
otherwise. 

Furthermore, if 7, J C {*, 1 , . . . , V + E} given as / = {ii < ■■■ < i k },J = {j\ < ■■■ < j k }, then we 
define the matrix L[7, J] to be L, with columns i\, . . . , i k removed, and with rows 71, ... , jk removed. 
It is straightforward to see that, if I = {*, 1, . . . ,E}, then L[I,I] = Q, where Q is as in eq. (14.631) . so 
because Uq = detQ, 

U G = detL[/,7] . (5.77) 

Obviously, since each matrix entry of L[I,I] is linear in the Schwinger parameters, it follows that 
Uq is a homogeneous polynomial of degree V. The polynomial 7g can similarly be written in 
terms of determinants of the Laplacian. Let us define, for i {*, 1, ...,E} 



Ii = {*,!,..., E,i} 



(5.78) 
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The matrix detQ is of course formed from (— 1)' +J times the ?'j-minors of Q, which are given 



by (—l) l+J detL[Ii,Ij]. So the polynomial Jq is the graph polynomial Jq is given by 



Tg 



det£ 



V+E 

£ (Q-^ijBi-Bj + C 

i,j=E+\ 



V+E 



£ (-l)' + MetL[/ i ,/ J ] £ aa&flZu 

i,j=E+l 



k,le{l £}: 

(*),C//)eEG 



+ detL[/,/]£ £ 



a 



i=l J€{1,. ..,£}: 

(.7) e EG 



(5.79) 



The determinants of L[I, J] are closely related to the trees within the graph G* by the so-called 
(generalized) "tree-matrix-theorem" ll56l 171142 1 . To state this theorem, let us first give some nota- 
tion. A tree T in G* is a connected subgraph of G* without loops. It is called a "spanning tree" 
if it has the same vertices as G*. A forest F is a collection of disjoint trees, i.e. a subgraph of G* 
containing no loops. We denote by the set of all forests, and by ^ the set of all forests 



{T u ...j k )e3T k 



(5.80) 



consisting of precisely k trees. Such a forest is called a "fc-forest". A forest F for the graph G* 
above is drawn in the following picture: 




•X m a forest FeT 4 ( 1,2) 



With a forest such as F, we associate a monomial mp formed from the product of all Schwinger 
parameters associated with each line. In the above example, 

m F (a) = Y[ a ki = «i5 oc 59 OC910 a 96 a*7 0C84 ocsio • (5.81) 

(kl)eF 

The graph polynomials Uq-, 7g are linear combinations of such monomials, where F runs through 
certain forests of G*, as we will explain. 

If I, J are as above, then denote by 3? k I,J the set of all k spanning forests in G* with the property 
that each tree T within a forest F G £? k I,J contains precisely one vertex i a from / and one vertex j'p 
from J. The forest clearly defines a bijective mapping from / — > J, i.e. a permutation, and we call 
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sgn n F the sign of this permutation. Finally, we call |/| = i\ + ... + i k , and similarly for 7. With this 
notation in place, the tree-matrix theorem states that 



detL[/,7] = (-l)l / l+l y l £ sign(7i F ) \[ a 



v 



(5.82) 



Applying this identity to the set / = 7 = {*, 1, ...,£} gives, with Te+i = ^^j, 

«c= £ m F (a)= £ EI a 0'- 

One can also obtain a formula of this nature for Jq- For this, we decompose 7g as 



(5.83) 



7g 



detg 



£ (Qr%B r Bj + C 

i,j=E+l 



(5.84) 



where the polynomial 1/q has been defined as Vq = ^q{{\}). We now apply the matrix tree theo- 
rem to each polynomial. Letting Xe(z, j) be the is-forests such that precisely one tree is connecting 
the vertices i, j, and such that the remaining E — 1 trees each contain precisely one of the points 
{*, 1, ...,£■} \ {/, j}, we can write 



? = £ (1-Z (7 ) £ H a«. 

l<&j<E FeZ E {iJ){kl)€F 

Furthermore, we obtain, taking into account various cancelations of negative terms, 

E 

^=£ £ n a kl- 

i=lFe1 E {i,*){kl)eF 

Putting these equations together, we get 

7 G = £ (l-Zij) £ Yl a kl + £ £ a kl 

l<i^j<E FeZ E {iJ){kl)eF l<i<E Fe1 E (i,*)(kl)eF 



(5.85) 



(5.86) 



(5.87) 



To illustrate the formula, we give the polynomials Uq-, 7g when G is the setting sun graph. The 
graph in question is: 




The corresponding graph G* is: 
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The external points are X\,X2, and the only point-pair invariant is hence Z\2- The only spanning 
forest in F e £2(1 , 2), drawn in red is: 



X\ X3 X4 X2 




* 



The other 8 spanning forests in F G I2 '■= ^seMM^li^s) are: 



X\ A*3 X4 Xi X\ X3 X4. X2 




X\ X-i, X4 Xi X\ X3 X4 X2 




X\ Xt, X4 Xi X\ X3 X4 X2 




X\ X3 A4 A*2 X3 A4 




There are also 6 spanning forests F e T3. The graph polynomials are, according to (I5.83I) . (I5.87I) : 

i/g = OC13OC3* + ai3a34 + ai3a42 + a42a2* + OC23OC2* + OC23OC3* (5.88) 

7 G = (1 -Zi 2 )ai 3 a 2 4a34 + 

+ ai3a3*a4* + a^as^a^ + ai30C34a4* + ai3(X3*a24 

+ a3*a4*a42 + 0634064*0x42 + 063*06340642 + 061306*40642 (5.89) 

Thus, in summary, we have shown in this section that: 
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Theorem 3. • The graph polynomial lie is homogeneous of degree V in the Feynman param- 
eters CLij,i,j G {*, 1, . . . ,V + E}. Each monomial has coefficient +1. The monomials are 
related to spanning forests of the graph G*, from the set 

%e+i — {spanning forests of G* with E + l trees, 

each containing precisely one vertex from *, 1, . . .,E} , (5.90) 

see eq. (15.831) . 

• The graph polynomial Jq is homogeneous of degree V + 1 in the Feynman parameters 
a i/>*» j : e {*) 1» •• •) V+E}. Each monomial has coefficient +1 or (1 — Zy) =H {Xi— XjY/2, 
which is > for points X, G S D , but which can become negative for points G dSp. The mono- 
mials are related to spanning forests of the graph G* from 

%e= IJ ^M, (5.91) 

r,s€{*,l,...,E} 

where Ze(t,s) is the set of all spanning forests with E trees such that each tree contains 
precisely one vertex from *, 1, . . .,2s, except for one tree connecting the vertices r,s. The 
formula for f is eq. (|5.87l) . 

For completeness, let us record the expression for fa in terms of other coordinate systems, 
which are obtained by simply substituting the corresponding expressions for the point-pair invari- 
ant, see table El Letting the external points X{,i = l,...,E be at equal time, i.e. = (Y,x ; ), the 
expression for Jg in the cosmological chart is: 

j G =H 2 z mt £ (x ; -x,) 2 £ n ««+ e n a « ( 5 - 92 ) 

\<i<j<E Fet E (i,j) {kl)eF Fe% E (kl)eF 

Letting the points X;, i = 1, . . . ,E be at equal time in the static chart, i.e. = (r|,r/,JC(), and letting 
Hri = cos 0,, 0, G (0, 7l) , the expression in the static chart is 

7g= £ (sm 8; sm fy + cos 6/ cos 6y- *,••.£/) £ \\ a«+ £ £ ]^[ a H 

\<i^j<E Fe% E (i,j){kl)eF l<i<E Fe% E {i,*)(kl)eF 

(5.93) 

In terms of the signed squared geodesic distances a, 7 , the polynomial is obtained by replacing 
1-Z ;7 = 2sin 2 (f v ^7). 

5.2 Mellin-Barnes representation 

Using the structure of the polynomials U , f derived in the previous subsection, one can derive 
another "Mellin-Barnes-type" representation of the Feynman integral Iq, which now describe in the 
present section. As in the previous sections, we will assume in this subsection that all points X( G S D 
are in the Euclidean section of the complexified deSitter space, meaning that 1 —Zn G (0, °°) for all 
i, j = l,...,E. The Mellin-Barnes representation will allow us to obtain the analytic continuation 
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to dSo, as we explain below in subsec. |5~4l and it will also allow us to derive the exponential decay 
of the correlators below in sec. [6J 

To state our result, we recall that the monomials in 11 resp. J were labeled by forests F from 
X£_|_i resp. %e- Each such forest is a union of E + 1 resp. E trees within a graph G* that is 
canonically associated with G, see the previous subsection for the precise definition. There is one 
distinguished trivial forest in this collection, called F = 4>, given by 



$ = {((£ + lH...,((£ + y»}6%i 



(5.94) 



Our Mellin-Barnes formula will involve an integral over multiple paths in the complex plane over 
variables wp. There is one such variable for each forest F, except for the forest F = 4>. Further- 
more, it involves the functions H n defined by the multiple integral 



H„(z) 



-2- 



'n-l 

n 

.1=1 



! '°° dsi 



(5.95) 



n-l 



n-l 



n-l 



n-l 



T(^ + ip+z-^ Sl )T( V -/p+z-E s t )r(-z + £ si)T(-z + £ *, - 



1=1 



1=1 



n-l 



n r(^i + ip + 5/)r(^i - /p + Sl )r(- Sl )r(- Sl 



i=i 

D-2\ 



1=1 



1=1 



The ^/-integrations in H n (z) are along paths parallel to the imaginary axis in the strip (IC.148I) . 
H n (z) is analytic in the strip < (R( z ) < -2^1. 

The somewhat lengthy proof of the following result is given in appendix O 

Theorem 4. The Feynman integral Ig(Xi, . ..,Xe) for non-coinciding Euclidean configurations of 
external points X\, . . . ,X E G S D is given by 



I G = K G f r G (w) J] (l-^) 1 ^, 

«/ W 1 £ „^ L 



l<r^<£ 



(5.96) 



where the sum £ F in the exponent is over all forests in %e(t,s), i.e. those connecting the vertices 
X r and X s , see the following picture for an example. Kq is the numerical constant in eq. (|4.71l) . and 
Z rs are the point-pair invariants H X r -X s . Furthermore, Tg(w) is the meromorphic kernel 



H^ + EfWf) Uf T(-W F ) Il(y)^4 ^(£F3(y)WF) 



r G (vv) 



'D+l 



(5.97) 



with riij > 1 the number of lines between a pair of vertices in G, and with the forest <J> omitted 
in any sum/product over forests F . The integration f~ = YIf^® I i s over paths which are 
asymptotically parallel imaginary axis which: 



• leave the poles of r(— Wf) < for F ^ 4> to the left, 



36 




series of poles of H njj (Y,FB{ij) w f) to the left and the +Nq 



Eq. (15.961) is one of the key results of this paper. It expresses the Feynman integral Iq as a 
multiple Mellin-integral. 



The Feynman integral Iq is of the form of a so-called "(generalized) //-function" (see appendix lAl). 
and the conditions for absolute convergence of the contour integrals over w stated there are fulfilled 
in Iq, as we also show in appendix|Cj This implies in particular, as we will discuss in more detail in 
subsec. l5~4l that Iq(X\ , . . . ,Xg) defines an analytic function (with cuts) in a subset of configurations 
in the complex deSitter spacetime. 

5.3 Example 

We will now illustrate the formula (15.961) for the deSitter space Feynman integral Iq in the case of 
the simple graph G with E external lines and V = 1 given by the following picture for E = 7: 

x, 



The integration variable is Xe+i, which is X& in the above example. For E = 4, the corre- 
sponding Feynman integral Iq gives the first order perturbative correction to the 4-point function 
(<K-Xi) • • • ( r , (^'4))o,?i- As always, the corresponding graph G* has one more vertex * and is shown 
in the following picture: 





x 1 T x 2 
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The possible forests in the sets %e an d T^+i except for $ = {(*(£'+ 1))} are depicted in the 
following pictures. There is one integration variable Wij,i,j £ {*, 1, . . . ,E} for each of the forests 
{(i(E + 1)), ((E+ connecting vertex i with j. Furthermore, there is one integration variable 
Wi,i G {1, . . . ,E} for each of the forests {(i(E + 1))} connecting vertex E + 1 with with i. The 
collection of these integration variables is called w as in the previous subsection. 



x 7 ' 


' x 2 


x 6 \ 


/ *3 








A forest in Xy(2, 6), corresponding to W26- 



A forest in %g, corresponding to w%. 



The kernel Tq in eq. (15.971 ) becomes 

^(^Y L +LkZk-Lk^iWki) Uk<i r {-m) Uk^(-wk)T(-zk + w k +Zi^kWki) 



r G (w) 



Yl 2 -2»r( V + 'P + a)r( V - *p + ^) r ( - ^¥ - *0 



(5.98) 



for this simple graph, where Zfe = Wjt + Wjt* + Y.i^k w ki- All summations/products are between 
1 , . . . , £. The resulting integral (15.961 ) can be simplified somewhat in this simple example, by chang- 
ing variables to Zk, and making repeated use of the formula^! 



dv r(b + v)T(d-\) _ r(c-b-d + e-\)T(b + d) 
2niT(c + \)r(e-v) ~ T(e-d)T(c + e-l)T(c-b) 



(5.99) 



valid for fR(c + e — b — d — 1) >0to carry out the Wk integrations, using also the duplication 
formula T(x)T(x + 1/2) = 2~ lx+l ^/iiT(2x). We do not give this somewhat lengthy calculation 
here but only record the final result, which is: 



Ig 



2 D K G 



^Uf r((D - 1)/2 + ip + zi)T((D - l)/2 - ip +zi)T(-zi -{D- 2)/2) 



k<l Jw kl 



T(D + Zizi) 



k<l 



i-z, 



■Id 



(5.100) 



In this expressions, all sums/products are again between l,...,E, and the z/-integration contours 
are parallel to the imaginary axis in the strip (IC 1481) . and the wm -integration contours are parallel 
to the imaginary axis with %K{yvu) < 0. The constant Kq is as defined in eq. (14.711) , with V = 1 
there. Our formula agrees with that found by 11451 using another method. 



19 



This is a version of the i /^-summation formula [cf. eq. ( IA. 1181 )1. written in Mellin-Barnes form. 
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5.4 Analytic continuation 

So far, we have focussed on the Feynman integrals Iq for non-coinciding (X; 7^ Xj) Euclidean con- 
figurations (Xi G S D ) of points, i.e. ones such that all point-pair invariants Z,j G [—1,1). To obtain 
the correlation functions (<f> (Xi ) . . .§(Xe))xo as distributions for deSitter configurations (Xj G dSp), 
we must, as already indicated above in eq. (|4.30l) . analytically continue through the complex de- 
Sitter configurations (Xi G dSp). We will outline this in the present subsection. 

The free field 2-point function (i.e. X = 0) was seen above in sec. 13.31 to be analytic on at 
least the domain T 2 = {(X h X 2 ) G (dS%) 2 | 3(Xi -X 2 ) G V + }. This domain by itself does not 
contain real deSitter configurations (Xi,X 2 ) G (dSr>) 2 ; those lie on the boundary of this domain. In 
fact, the free 2-point function is defined as the distributional boundary value as we approach this 
boundary, see item 2) of sec 13.31 and footnote [TT] The full domain of analyticity of the free field 
2-point function is larger than T 2 , and can be obtained e.g. by acting on T 2 with complex deSitter 
isometries from 0(D, l) c . Alternatively, we may observe that the two-point function is, where it 
is analytic, a function only of Z. As (X\,X 2 ) varies through T 2 , Z varies through the cut complex 
plane C\ [1,°°). The two-point function is analytic for any configuration such that Z is in this cut 
complex plane, and this includes many real deSitter configurations, such as when (X\,X 2 ) G (dSo) 2 
and spacelike. One sees more generally that, using Wick's theorem, the free-field it-point functions 
are analytic at least in the domain 

T E := {(X h ...,X E ) G (dS^) E I 3(Xj-X j+ i) G V+} . (5.101) 

Bros and Moschella (see ref. 1 of [9]) suggested that this should be contained in the domain 
of analyticity also for interacting quantum field theories. Again, real deSitter configurations 
(X\, . . . ,Xe) G (dSo) E are only on the boundary of 7g, and the it-point correlation functions on 
real deSitter are defined as distributional boundary values [9|. They also showed how the domain 
Te can then be enlarged using complex deSitter isometries, the "edge-of-the-wedge-theorem" etc., 
after which it also contains many real configurations. We will now consider whether T E is indeed 
a domain of analyticity for the correlators in the interacting theory, in the context of perturbation 
theory. Since the is -point functions (§(Xi) . . .§(Xe))q % are given in terms of the Feynman integrals 
Iq, we must analyze those. 

For a configuration of E points (X\,. ..,Xe) G Te, the point pair invariants are all in the cut 
complex plane 

Z„eC\[l,oo) for r^s. (5.102) 
Indeed, letting e.g. r < s, we can write 

s-l 

3(X r -X s ) = T3(Xj-X j+1 ) G V + , (5.103) 

Pr" v ' 

and by looking at the cases Ui(X r —X s ) G V+UV -- ,^ V + UV~ separately, we get from this 
H 2 (X r —X s ) 2 ^ [— 1, —00), so Z rs = H 2 X r -X s [I, 00 ). Even though there are many dependencies 
among the Z rs , the Feynman integrals Iq can be viewed as functions of the Z rsi r ^ s, considered as 
independent variables, which are defined when Z rs G [—1, 1). We must ask whether it can be ana- 
lytically extended to Z rs in the cut complex plane, for then Iq, viewed as depending on X\ , . . . ,Xe, 
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is also defined as an analytic function on Tp. For this, we can look at the explicit form of Iq as 
functions of the Z rs , viewed as independent variables. The most useful representation of Iq for this 
purpose is (I5.96I ). 

That formula gives Iq as a //-function, whose analyticity domain is described in appendix lAl 
According to eq. (IA. 1241) there, we must calculate the real numbers Ap associated with each inte- 
gration variable wp in eq. (15.961) . in terms of which the domain of analyticity is then (at least) 

|arg(l -Z rs )\ < ^%mf{A F | F e Z E (r,s)} . (5.104) 

Looking at the explicit form of the kernel Yq in eq. (15.971 ), we derive in appendix O step 7), that 
Ap > 2 for any forest F, from which it then follows that Iq is indeed analytic in the Z„ in the 
cut complex plane. When viewed as a function Ig(X\, . . . ,Xp) of the external complex deSitter 
configuration, Iq is hence analytic at least on Tp C (dSp) E , but even for all complex deSitter 
configurations such that the Z rs are all in the cut complex plane. This set contains many real 
deSitter configurations in (dS^) E , e.g. if all points are mutually spacelike. To define Iq on the 
entire space of real deSitter configurations, it must be understood as a distributional boundary 
value. 

In order to define this distributional boundary value, one can again make use of the represen- 
tation (|5.96l) together with a general result in distribution theory, see e.g. ch. IX of [|34l . This 
general result states that if u(z) is an analytic function that is defined on an open subset of C" of 
the form U + iC, where U C W 1 is open and C is an open convex cone in W 1 , and if u(x + iy) does 
not increase too fast as C 3 y — > in the sense that 

\u(x+iy)\ < est. \y\~ N for all y G C, someN, (5.105) 

then u possess a distributional boundary value for C 3 y — > 0. The same applies if u is only defined 
and analytic for y in a neighborhood of the tip in C. When u is defined on a complex manifold, 
then the boundary value is defined in each chart and then pieced together by a partition of unity. 
Now let {X\ + iYi, . . . ,Xe + iYp) G ( dS^) E fl Tp be a complex deSitter configuration, with X,-, Y,- 
real. Then, Y r — Y S E V + for r < s, and because 

|^ 2 |>|(3^) 2 | if^GV+, (5.106) 

we get |1 — Z rs \ > j\(Y r — Y s ) 2 \. Because Iq is a generalized //-function that is defined by an 
absolutely convergent multiple contour integral, we can use this to estimate for (X± + iY\ , . . . : Xp + 
iYp) in a bounded subset of T E : 

\I G (X 1 + iY h ...,X E + iY E )\ < K G J] |1-Z„| M *>^) /|r G (w)| (5.107) 

< est. J] \(Yr-Y s ) 2 \- N , 

\<r<s<E 

where the Yq is the integral kernel (15.97k and the sum is over all forests F in X£(r, s) connecting 
vertex r with vertex s. In the second line we have used the fact that the multiple contour integral is 
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absolutely convergent, and we have set —N = inf^ f $H(wf). Let Ce-\ C (M d+1 ) £ 1 be a convex 
cone with small opening angle around the vector (e, . . . , e) with e = (1, 0, . . . , 0) G M Z)+1 . Then, for 
(X 1 -Y 2 ,...jE-i-YE)eC E -uv/e have 

J] \(Y r -Y s ) 2 \- N <cst[\Y l -Y2\ 2 + --- + \YE-i-Y E \ 2 r M (5.108) 

l<r<s<E 

for some M, where | . | denotes the Euclidean norm of a vector. Combining this inequality with 
eq. (15.1071) then evidently gives 

\I G {X X + iY h . . . ,X E + iY E ) | < est. -Y 2 \ 2 + --- + \Y E - l - Y E \ 2 ]- M (5.109) 

for a bounded set of (Xi +iY h ...,X E + iY E ) G (dS%) E such that (Yi Ji - Y 2 , . . . ,Y E -\ - Y E ) G 
M. D+l x Ce-\- Expressing this result in any coordinate system x + iy of (dSp) E centered around 
(X\, . . . ,Xe) then leads to the desired inequality (15.1051 ) for a suitable cone C related to MP +l x 
Ce-i- 

Thus, the upshot of our discussion is that, for real deSitter configurations, Iq is the distributional 
boundary value of an analytic function. The boundary value prescription given above can be 
summarized as saying that [compare eq. (15.961 )1 

, c{ x u ...,x E) =2K [r c ^ n ( (J^±piM. f WF (5 , 10 ) 

J ™ l<r<s<E \ 1 J 

where the Xi are now in the real deSitter space, and where e = (1 , 0, . . . , 0). As usual we mean 
here that Iq has to be smeared with a testfunction f(X\, . . . ,Xe) first for £ > 0, and then we send 
£ — )• 0. The sum over F in the exponent runs as usual over all is-forests in the graph G* such that 
the external leg associated with X r is connect to that with X s by a tree in the forest. 



6 Quantum no-hair/exponential clustering 

We will now argue that the connected correlation functions have exponential decay in timelrke 
directions, and a corresponding decay in spacelike directions. As above, we consider a princi- 
pal series scalar field (c = — (D — l)/2 + ip, p G M), and we restrict attention to configurations 
(X\, . . . ,Xe) in the real deSitter spacetime such that no point is on each other's lightcone. For 
general configurations, we would have to take into the distributional nature of the correlation func- 
tions. This case could also be dealt with using our methods, but for simplicity we will leave it 
aside. 

Thus, we take one point, X r , and move it away from the other points to infinity either in timelike 
directions so that Z rs — > +°° for any s ^ r, or in spacelike directions so that Z rs — > — °° for any 
s 7^ r@ We will argue that the correlation functions satisfy, to all orders in the coupling constant 
I: 

(<S ? (X l )---$(X E ))l = O(\Z rs \-^ as any \Z rs \ -> °°, (6.111) 

20 Note that, since deSitter space has horizons, we could also let X r go to in such a way that X r does not become 
eventually time-like to any of the X s for s ^ r, and so in general it need not be true that Z„ — > °° for all s. 
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for any ^ < (D — l)/2 and any s = 1, . . . ,r, . . . ,E. Here, Z„ is the point pair invariant, which 
for timelike related points is given by proper time separation % rs between the two points X ri X s as 
Z rs = cosh//x„ ~ e^l 1 ™!, and similarly for spacelike related points, see eq. (12. 1 2b . For large spatial 
separation, e.g. in the cosmological chart, see (12.91 ), this implies: 

(§(t,xi)---§(t,x E )} < i. = O(\x r -x s \- 2 ^) asany |x,.-x,| -> (6.112) 

For a scalar field in the complementary series, our decay result (16.1 1 II) would hold for £, with ^ < 
(D — 1)/2 — y/(D— l) 2 /4 — m 2 /H 2 , and we would get a correspondingly weaker decay. Again, to 
keep things as simple as possible, we will not consider this case here. 

We demonstrate (16.1 1 II ) in the remainder of this section. Evidently, since the connected cor- 
relation function is a sum of contributions from connected Feynman graphs G, it is sufficient to 
demonstrate the exponential decay for each such contribution, i.e. 

I G = 0{\Z rs \-S) , (6.113) 

for \Z rs \ — > °o. Our main tool is the Mellin-Barnes formula (15.961) for Iq. The desired exponential 
decay (16.1 131) will be obtained by looking at the location of the contours in this formula. Indeed the 
Mellin-Barnes formula immediately gives for X r going to infinity in time- or spacelike directions: 

|/ G (Xl,...,X G )| <cst. J] |Z r ,| su P^^) L\T G (w)\ <cst. n |Z r ,| su P^^) , (6.114) 

where the sum is over all forests F G %E(r,s) in the graph G* consisting of E trees, one of which 
connects the vertex X s with X r , see subsec. 15.11 for the precise definitions. The contour integral 
is absolutely convergent, as discussed in sec. (|5.4I) . The supremum is taken along the integration 
paths followed by the variables w, see subsec. I5T21 These contours can be somewhat complicated, 
so to get an estimate of the supremum, it is best deform the contours followed by the wf for such 
forests to straight lines parallel to the imaginary axis. There are many ways of doing this. 

Let us first suppose, for ease of notation, that E = 2 and r = 1, so that s = 2, and let us also 
suppose, that a pair of vertices in G is connected by at most one line [so that n ;j = 1 for all i, j in 
eq. (15 .971) 1 . Take any forest i , GT 2 (l,2), J Pn4> = [cf. eq. (15 .941) 1 and move its path of integration 
to be at 9\(wp) = est. with 

~^-<D\(w P )<-^-+£ . (6.115) 

For each (kl) P, (kl) 4>, choose a forest F(kl) G T 3 containing (kl), but F(kl) HP = 0. We 
deform the integration path for the corresponding variable Wpm\ so that it is the same as that for 
wp. We deform all other wf's so as to run along paths with — £ < %K{wp) < 0. The paths have been 
chosen such that any pole at No of T(— wf) remains to the right, such that the — No±/p — (D— 1)/2 
series of poles of each H\ (Ef3(W) w f) remain to the left and the +No — (D — 2)/2 series of poles 
remains to the right of the integration paths. Furthermore, sup{9^(w/7) I F G ^(l^)} < — £, for 
sufficiently small thus showing exponential decay. 

However, when deforming the contours, we will cross some of the poles of T(^±± + £ F w F ) 
to the left. To avoid this from happening, we choose another forest Q G X3 and we translate the 
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contour followed by wq to the right by the appropriate amount. The price that we pay is that now 
we cross some poles of F(—Wq), and we pick up corresponding residue. Each such pole term 
looks exactly like (15.961) , but it has no integration over wq, which is set to some n E Nq. There is 
one such term for each pole crossed. But since all the variables Wf appearing in the exponent of 
the factors (1 — Z rs ) are associated with F's from %i, they are not affected by this, as Q E As a 
consequence, the pole terms can be estimated in the same way as before. The general case is dealt 
with in exactly the same way; without loss of generality we may again assume that r = 1,5 = 2. 
The forests F(kl) 1 Q are now in 1e+i- When some ntj > 1, a similar argument can be made. 



7 Conclusions and outlook 

In this paper, we have: 

1. Given parametric formulae for Feynman integrals for massive scalar fields in deSitter space- 
time for an arbitrary graph G, see eq. (15.961) . These integrals are the building blocks of the 
perturbative expansion of deSitter correlation functions (in the Bunch Davis state = Har- 
tle Hawking state = Euclidean vacuum state). The parametric integrals involve a multiple 
contour integral over parameters which are in correspondence with tree graphs, and more 
generally forests, within the graph G. An alternative parametric representation in terms 
of certain graph polynomials associated with G was also given in eq. (14.691) . Our para- 
metric representations have certain features in common with similar expression known in 
Minkowski quantum field theory. However, they are different at least in that (a) we work in 
position space and (b) our forests and effective graphs associated with G involve a "virtual 
vertex", *, that is, at some level, a reflection of the fact that the loop integrals are restricted 
to the deSitter hyperboloid. We were not able to find a "momentum space version" of our 
formulae, mainly due to the fact that no useful and simple momentum space expression for 
the Feynman propagator appears to be known in deSitter. 

2. We have found that the Feynman integral belong to a class of special functions called "gen- 
eralized //-functions" ||49l , which are a generalization of the hypergeometric function. 

3. We have discussed the renormalization and convergence of our parametric form of the Feyn- 
man integrals. 

4. We have discussed the analytic properties of the Feynman integrals/correlators, in particular 
how to go from Euclidean to Lorentizan configurations via analytic continuation. 

5. We have used the parametric representation to show the exponential decay in timelike di- 
rections of the deSitter correlation functions, to arbitrary orders in perturbation theory. As 
we have explained, these results show that the correlation functions of essentially any state 
will approach that of the vacuum at late times, at an exponential rate. This statement may 
be viewed as a quantum analog of the cosmic no hair theorem, for a test scalar field with a 
positive mass. 
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It would be very interesting to try to generalize our results in the following directions: 

1. Look at exponential decay for other types of fields, in particular to massless fields such as 
(i) massless scalar/spin- ^-fields (ii) Yang-Mills-fields (iii) the graviton. For the free spin- 
^-propagators, see |fT6ll , for a general discussion of the renormalization theory of the Yang- 
Mills-field on general curved spacetimes see 11331 . for free Yang-Mills- and ghost- in deSitter 
see If28ll . and for free graviton propagators see e.g. 11551 l26ll28Tl . 

2. Look at the decay of correlation functions in the Schwarzschild-deSitter spacetime, depicted 
in the following conformal diagram: 




J»-3 

d s 2 = -f dr\ 2 + f- 1 dr 2 + r 2 d(Q 2 D _ 2 , f=l-^—-H 2 r 2 . (7.116) 

In the gray region, one would expect the decay of the fields/correlators to be essentially the 
same as in pure deSitter spacetime. But in the static region, the influence of the both the 
black hole and cosmological horizon will be felt, and the likely outcome is less clear. A de- 
cay of the correlation functions in time-like directions (|r|j — T|j| — > °°) in this region would 
indicate a kind of quantum stability of the black hole, at least in the test-field approxima- 
tion. Unfortunately, a systematic study along the lines of methods of this paper seems not 
possible, because even the propagators of the free field theory in this background are not 
known in sufficiently explicit form. However, it might be possible to extract enough infor- 
mation about their asymptotic behavior e.g. from known properties of quasi-normal modes 
in Schwarzschild-deSitter IfTTl . 
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A Hypergeometric functions and their generalizations 



In this section we give formulas for the Gauss hypergeometric function and the free propagator 
referred to in the main text. The hypergeometric function is defined as 



2 Fi(a,b;c;z) 



v (a)n(b) n 

L — z 



n>0 



n\{c), 



(A. 117) 



where (a) n = T(a +n) /T(a). The series converges and hence defines an analytic function in the 
open unit disk in the complex z-plane. A value used in this paper is 



2Fi(a,b',c,-l) 



T(c-a-b)T{a)T(b) 



T(c-a)T(c-b) 
A transformation formula used in this paper is, for m = 1,2,... 



(A. 11 8) 



, , s T(m)T(a + b — m) , , m %-} (a — m) n (b — m) n 

/ l)m r(a + b-m) y (a) n (b) n 
1 ' T(a-m)T(b-m) ~ Q n\(n + m)\ 

x[log(l -z) -\\r(n + l) -\\r(n + m+l)+\\f(a + n)+\\f(b + n)](l -zf (A.119) 

where \|/(z) = r'(z) /r(z). Application of this formula to the vacuum 2-point function of the free 
Klein-Gordon field (|3.26l) gives for even dimensions D: 



D-2 



(4tc) d / 2 ) V 

x D/ y 2 (-c-D/2+l) n (c + D/2) n 
X to n\(2-D/2) n 



H 1 {X l -X 2 f 



-D/2+l 



(A. 120) 



H 2 (X l -X 2 ) 



2in 



+ 



+ 



1) 



D/2 



T(-c + n)T(c + D-\+n) 



r(-c-D/2+l)r(c + D/2) 



n>0 



nT(D/2 + n) 



H 2 (X { -X 2 f 



x [log — { — A — y~ - + 1) - \|/(n + D/2) + V|/(-c + n) + V|/(c + D - 1 + »)] 



The following asymptotic formula is known for large values |L| when |arg(Z— 1)| < n, |arg(L)| < 
71, see 11381 : 

x{ I c _ l ((L+(a-b)/2)%) + 

+I C - 2 ((L +(a-b) /2)x) [(c - 1 /2) (c - 3/2) (x^ 1 - cotx)] /{2L + a-b) + 
+(2c-a-b-l)(a + b-l)tanh(x/2) } (A.121) 
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where coshx = Z, and I v are the modified Bessel functions. This asymptotic expansion is valid 
uniformly in Z in the indicated domain. 

There are many generalizations of the hypergeometric function. A very general class is pro- 
vided by the so-called (generalized) //-function of several variables, which is defined by the 
Mellin-Barnes-type formula 

f (A) (B) (a) (b) 7 \-(u ( dwk 

H \{C) (D) (c) (d)> Zh ---> Zn )-[l\J Kk 2ni 



UU r( aj + LtAijWi) n y ii r( i-bj- LiB 



tjWi) 



Uj^ncj+ZiQjWi) n5=ir(i - dj-ZiDijWj, 



zp--<". (A. 122) 



The paths of integration Kj are intended, if necessary, in such a way that all the poles of T(aj + 
YsiMjWi) are separated from the poles of T(l — bj — Y^iBijWi). The parameters A, ; -,5,j,C,- 7 -, A; G K 
are all assumed to be > 0, and a^b^c^di G C. A discussion of such kinds of functions including 
functional identities, asymptotics, and their relation to various other classes of special functions 
can be found e.g. in [49]. The above integrals converge absolutely and define an analytic function 
(at least) when the n are not equal to zero and when 

it. A- 

|arg( Z; )|<^ (A.123) 

where 

V=£Ay-£Q+£2?y-£Dy . (A.124) 

.;' j j .i 

The formula reduces to the Gauss hypergeometric function and standard simple generalizations 
thereof in special cases. 



B Spherical harmonics, Gegenbauer polynomials, Kallen- 
Lehmann representation 

Here we give a "spectral representation" of the 2-point correlation function of an interacting field 
analogous to the Kallen-Lehmann representation in Minkowski spacetime. The derivation of the 
formula involves spherical harmonics in D-dimensions, so we briefly recall their basic properties, 
see e.g. J58H21 for more details. 

Spherical harmonics on the unit S D can be introduced via harmonic polynomials in the em- 
bedding space M D+1 . A polynomial P(X) on R D+l is called homogeneous of degree h if P(XX) = 
A, P(X), and it is called harmonic if it is a solution to the Laplace equation on R. D+l . The harmonic 
polynomials of degree h = L form a vector space of dimension N(L,D) — ( 2L +^^)(^+^ 2 ) ! - spher- 
ical harmonics on S D of order L are by definition just the restriction of the harmonic polynomials 
to S D . The spherical harmonics Yu(X),j = 1, ...,N(D,L) may be normalized so that 

Y, Y Lm(XiYY Lm (X 2 ) = 8(Zi,Z 2 ) , f dfi(X) Y Lm {X)*Y L , m ,(X) = 8 L ,A, m > (B.125) 



46 



where the 8 function is that on S D , defined with respect to the measure d/u. Expressing the Lapla- 
cian on M° +l in polar coordinates, on sees that the spherical harmonics are eigenvalues of the 
Laplacian V 2 on the D- sphere with eigenvalue —L(L + D— 1), so that L may be viewed as the ana- 
log of the total angular momentum, and m may be viewed as the analog of the magnetic quantum 
numbers. One has 

N(D,L) 9r_i_n_i 
I YUXiTYUXi) = Zu, D .n C [ L D - l)/2 (Z) , (B.126) 

where are the Gegenbauer polynomials, where Z is the point pair invariant. The Gegenbauer 
polynomials are expressible in terms of a hypergeometric function, 

c (D-i)/2 {)= r(L+D-i) (_ LL+D _ UD/2 A_^l) (B127) 

L V ; T(D)T(L+1) V 1 2 ) V 

Eq. (|B. 1 261) may be viewed as saying that the Gegenbauer polynomials are, up to normalization, 
the integral kernels of the projector onto the eigenspace for the eigenvalue —L(L + D— 1) of the 
Laplacian on S D . Since the dimension of this eigenspace is equal to N(D,L), one gets the orthog- 
onality relation 

f dZ (1 -Z 2 )^ 2 - 1 C^ 1)/2 (Z)C[^ 1)/2 (Z) =N DiL h.L> (B.128) 

for L,L' E No, with normalization factor 

vo\(S D ~ l ) 

Nd,l= vol(5g) N(D,L) (2L + D-1)- 2 . (B.129) 
By the same argument, one gets the formula 

- ^lyl^"^ -^^-^^-^!) <B ' 130) 

for the Euclidean Green's function of (— V 2 +m 2 ) on the sphere S D of radius H , where c is as 
in the main text. The above sum can be converted to a contour integral over L with the help of a 
Watson-Sommerfeld transformation, as observed in Ii44l : 

(<KXi)<KX 2 )) = f ^-.(2L + D-l)P L A L (Z) (B.131) 
Jc 2ki 

where the contour C is running parallel to the imaginary axis, leaving the poles in the denominator 
of 

% 1 / 2 T(D/2) 1 

Pt [ ~ T(D/2 + 1/2) -c(c + D- 1) +L(L-D+1) (B ' 132) 
to the left, and the poles at No of to the right. The kernel 

, s H D 2 T(L + D-l)T(-L) ( , 1+Z\ 

MZ) = m o /2 r(D/2 V 2F, ^—L,L + D— l;Z>/2; - J (B.133) 
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is equal to the free Euclidean Green's function for the mass parameter M 2 = —H 2 L(L+D — 1). Pi 
is interpreted as the "power spectrum", or "spectral density". The original sum is recovered if the 
contour integral is evaluated by means of the residue theorem; convergence of the contour integral 
follows from eq. (IA.121I) . Such representations are closely related to the so-called Kallen-Lehmann 
representation in Minkowski-space, and have been pioneered by [0 in deSitter spacetime. 

A similar representation is also possible for the two-point function of the interacting theory. 
In the context of perturbation theory, the analog of the above spectral representation (IB. 1 33b for 
the 2-point function (§(Xi)§(X2))qx for an interacting field has spectral density given by a sum of 
contributions from individual Feynman diagrams G, 



G 



= H "^77 sym(G) Pl.g , (B.134) 



i.e. Pi g is the contribution from an individual Feynman diagram G. A formula for this can 
be obtained straightforwardly using the orthogonality of the Gegenbauer polynomials, and our 
parametric formula (15.961 ) for Iq. To state the formula, let us define the graph G to be the graph 
obtained from G by closing off the two external legs, so that G will have no external lines, and so 
that it has one more internal line, called Let Ig(L) be the expression (15.961) , with one c — > L in 
the i7„ 12 -factor corresponding to the line I which we closed off. Then we have 

n l ' 2 Y{D/2) sirnzL (2L + D- 1) 7 G (L) 
L ' G L '° T(D/2+l/2) % -c(c + D-l)+L(L-D+l) { ; 

The contour C in the spectral representation must now run parallel to the imaginary L-axis asymp- 
totically, and it must separate the poles of Pi \ from those of Ai (at No). The spectral representation 
provides an alternative way of performing the analytic continuation of the two-point function from 
the sphere to deSitter spacetime. 



C Proof of thm. i 

We here give the lengthy proof of thm. |4] Recall from subsec. [5T| that 11 , J are polynomials of 
parameters a (J , where there is one parameter for each edge (ij) e G*, with G* the graph associated 
with G. In this appendix, we find it convenient to set 

a,= / e "' iO' = V = £ + l,...,£ + V 
I e V/ > otherwise if i, j ^ *, 

and we also set H = 1 . The parameters v ;j G R are real. The main tool in the proof will be the 
well-known Mellin-Barnes (MB) identity 



( Al + ...+A iV+1 ) (B (C.137) 
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for 9t(co) < 0, A,- > and for N contours Ki going parallel to the imaginary axis such that 

tt(w,-)<0, JK[a>-£wj] <0. (C.138) 

After these preliminaries, we come to the proof, which is divided into several steps. 

Stepl: For^(co) < -V(D + l)/2,0>«R(zy) > -(D+l)/2 the following integrals are absolutely 
convergent: 



L 



1 CO 



71 



r(-to) w(w,v) a) +( v + 1 )( Z) + 1 )/ 2 



r(-cp+Ejz;) n 7 r(-z, ; 
rf-co) 



(C.139) 



^ v nxj-x k y 

j,k 



(O-ZjZj 



Yl e " jZj ( x j 



2 Yi 



Here, there is one complex variable Zj for each variable uj, and the integral j x denotes an integral 
over (]R D+1 ) 1/ against the variables Xe+i, . . . ,Xe+v- The second equality sign, follows from the 
MB-formula, whereas the first one is obtained e.g. by manipulations similar to those in sec. 14.31 
Because the integrals are absolutely convergent, we may change the order of the integrals in the 
last expression. Multiplying also both sides by r(— co), we trivially get 



7I ^+i) r( _ Q) . 



V{D+1). f(u,v 



,co+V(D+l)/2 



(C.140) 



^(Xj-xty 



1 



Step 2: We now apply an inverse Fourier transform / due mx to both sides. The integrals are 
absolutely convergent and given an analytic function in x for 



D+l 



> 3(xj) > 



(C.141) 



and we get, redefining also CO — £ y ixj — > co: 



V(D+l) 



TZ 



V(£»+l) r ( 03 ^j ix J 2 



r(-co)n i r(-a J ) 



? (u,v) w+L j ix j+ v ( D + 1 )/ 2 

u ( U ^0)+Zjixj + (V+l)(D+l)/2 



(C.142) 



2 \ ix i 



This equation is valid now for 



^(Xj-X k ) : 



*R(co) < 



(C.143) 
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and for 3(xj) in the above range. 

Step 3: We integrate both sides of the last equation against (2ti)~ v f dx over a contour of constant 
3{xj) in the above range. We also use that for such a contour / dxj (Xj)' x j = 2% 8(Xj — 1) in the 
sense of distributions. We are allowed to use this distributional identity because the integrand on 
the right side is in L 1 UC°° w.r.t. the variables Xe+i,. ■ ■ ,Xe +v . Using also that the measure on the 
sphere S D is dfi{Xj) = 2~ l 8(Z? - 1) d D+l Xj, we therefore get 




rT( ca Y ix v{D+l 
D+i) / M 10 Ljixj 2 



CO 



(C.144) 



7 (u,v) a+ Zj ix j+ v ( D + 1 )/ 2 



r(-(D)n^r(-a y ) Jn (a ^^+Ljixj+(v+i)(D+i)/2 ' 

which is valid at least for <K(co) < and < 3(xj) < (D + l)/2. 

Step 4: For < it is legal to use the MB formula on the integrand on the left side, which 



gives 



^(Xj-X k ) 

J,k 



CO 



Uj, k n-z jk ) 

r(-to) 



tl (X x 



,x 



V+E 



(C.145) 



Here, there is one integration variable Zjk for each v; k , except for one distinguished but arbitrarily 
chose edge (jk) = e E EG*. We must have £ft(z/jfc) < along the integration contours for all zm, 
and Y.kjZki = G>. t% is as in eq. (|4.34l) . We now take Y[j J dfi(Xj) of both sides of the equation. 
For sufficiently small £ > and — £ < fR(zjk) < 0, the left and right side are in L 1 with respect to 
the integration over Xe+i,. ■ ■ ,Xy + E e S D . The integrals over z on the left side is also absolutely 
convergent due to standard estimates on the T-function (Stirling's formula). Thus, we can exchange 
the integrals over X and z, and we obtain 



-I CO 



J* 



(Xj-X k 



Uj, k n-z jk ) 

r(-to) 



Y[] sD dii{Xj) ) t t (X h ...,X v+E ) 



We take an inverse Fourier transform j dv e lv ' y of this, where there is one variable y # for each of 
the variables vjk- Then using the result of the previous step 3), we get 



Ylf^dfiiXjU t#{X h . . . ,X V+E ) = 



(C.146) 



2 -v n v(DA 



V(D+1) 



n r r r (- Lj ixj - Lj,k iyjk - -4 



_ _ <r ( U $)Ljixj+L jtk iyj k +V(D+l)/2 
-iu-x—iv-y J \ i ) 



hJx 1 1 7 T( avlH^r! iy jk ) k~ ufr^J+ZjMyjkHV+W+W 
which is valid at least for £ > 3(yj k ) > 0, > 3(xj) > 0, and for 

V(D + 1) 

;j{xj) > 

j 



^3(y jk )+^3(xj)> 

j,k 



(C.147) 
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Step 5: We now change the order of integration J~ L — > L L. This will be justified a posteriori 
in the end, because all integrals in the final expression will turn out to be absolutely convergent]. 
We then get the Fourier transform j U jQ-™-l-iv-y f a+v{D+\)/2 / u u+{v+i){D+i)/2^ where now ro = 

E^/'/t + E We next rewrite this by bringing the integrand of this into the form of a Mellin- 
integral, using the MB-identity, so that the Fourier transform can be read off. The MB-identity can 
be applied provided that exponents of both 11 , J have negative real part, i.e. that 

<*(»)< (C ,48) 

which is compatible with our previous condition (IC.147I) . There is basically one Mellin-Barnes 
integration parameter for each term in 11 , J . As explained above in subsec. 15.11 the monomials 
from J are labeled by elements F of the set 

%e= (J $e(U)- (C149) 
!'je{*,l, ...,£} 

Thus, F E %e if it is a spanning forest with E trees. One of these trees connects a pair i, j E 
{*, 1, . . . ,E}, whereas the remaining E — 1 trees each only have precisely one element from this 
set. There is one independent integration parameter wp for each such forest F, except for one 
dependent parameter corresponding to one particular forest. Similarly, the monomials in 11 are 
labeled by forests F E %e+u see e£ l- (15.831) . There is again one Mellin-Barnes integration parameter 
for each such monomial, except for one forest, which is dependent. The two dependent parameters 
are related to the independent ones by the formula 

ItVF = X(^±2) +m , (c.,50) 

for $ and 

£ „ F = - C+'f+'> - m , (C ,51) 

for 11. Furthermore, for any forest from either %e or %e+u we must have D\(we) < 0. Let us 
define the quantities 

1rs(w) = W F ' (C.152) 



and 



ixj(w) = ^ Wf , (C.153) 

F: F3(*j) 

iyjkW = L w f- < C154 ) 

F: F3(jk) 



21 Indeed, we may replace the numerator in eq. ( lC.146t by r(— i(l +B)Xj — J^j k^jk — V 2 ^ s ma k e 
all integrals absolutely convergent by eq. JA. 124b , so we are free to exchange the order of integration. In the end, we 
get absolutely convergent integrals depending on e. At this stage, we can take e — > under those integrals, the limit 
being uniform. We omit the details. 
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Then using the MB -identity, we get 



«(«,v) ro +( ,/ + 1 )( D + 1 )/ 2 
= L - ^Zol) , fr 1 -Z w )^exp[g-g+,y • v] ■ (C.155) 

■ / wr(CD+ i % L ) r +s 

To read off the Fourier transform, we need to convert the integral over w into an integral over the 
quantities x,y as defined in eq. (|C.152I) . and a set of complementary variables which we call h. 
Thus, we are looking for a map (x,y,h) i-> w(x,y,h) which must satisfy the following properties: 

1. From eq. (IC 1 531) together with eq. (1C.150I) . we must have 

£ w F (x,y,h)=V(D + l)/2 + Y,ix j + 1 £iyjk- (C156) 

Fe1 E j j,k 

2. From eq. (IC 1 531) together with eq. (IC 1511) . we must have 

£ WF (^h) = -(V + l)(D+l)/2-Y j ix j - y £iy jk . (C.157) 

FeZ E+ i j j,k 

3. From eq. (|C.153I) . we must have Xj = Lf9(*j) wf(*,J,«) and y^ = T 1 F3(jk)^F(x,y,h) and 

4. The map (x,y,h) i— > w(x,y,ti) should be linear in h, and have maximum possible rank com- 
patible with the previous constraints, i.e. \%e U ^e+i \ — |EG*| — 1. 

There are many ways choosing such a map. One way is to pick, for any (ij) e EG*, a suitable 
corresponding forest <J> ;; e T^+i such that 4> i; - 9 (zj), and e.g. one forest Fq G %e- Then the 
variables h are defined as those variables w F such that i 7 ^ {Fo,4> ;j - 1 (zj) G EG*}, and the variables 
x,y are as defined in eq. (IC 1 521) . The choice of the forests must be made in such a way that 
the transformation (x,y,h) — > w is invertible (it is automatically linear). It is not difficult to see 
that there are indeed (many) ways of achieving this. Furthermore, for a suitable choice of the h 
contours, the x and y contours may be chosen so as to be at £ > 3(yjk) > 0, > 3(xj) > 0, and 
this is compatible with the range for which eq. (|C 1461) is valid. Hence, we can read off the Fourier 
transform as: 

r (t (ti rf\®+V(D+l)/2 

r(co+^) [ Q -^-y^n — L__ 



r n F n-WF(x^h)) n m%) 
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Then it follows from eq. (IC.1461 ) that: 



YlJ^XjU t if (X h ...,X v+E ) =2- v n v ^J_ i 



Ujr(-ixj)Uj,k r (-iyjk) 



x 



This formula is valid for e > 3(y jk ) >0,^±> 3(xj) > 0. 

Step 6: We next analytically continue this result to a larger set of y. From sec. 14.21 we know 
that the left side is analytic in iy for which A^(iy) H —No = ©• When we analytically continue 
the right side, we must move the x, h contours to new contours in such a way that (i) the poles of 
r(— w F (x,y y h), (ii) the poles of r(— ixj), and (iii) the poles of T{iyj k ) all stay to the right. After 
that, we can assume that we have continued to y in the range 

n jk (D-\) , . n jk (D-2) 

— - > tyjk) > J \ — " • (C160) 

We multiply the result of step 5) by Y\j, k H njk {iyj k ) [cf. (|5.95l) l and integrate over all yj k , along 
contours in the above range. The result is 

h T(Zj ixj + r M iy it + ) y 

The left side, is up to a constant, the Feynman integral Iq, by eq. (14.331) . 

Step 7: In order to get the desired formula for Iq, we need to change the variables of integration 
from x,y,h back to w, which by definition consists of all wp,F G %e UT^+i, except for the forest 
F — 4> given by 

<& = {((E + l)*),...,((E + V)*)}eX £+ i (C.162) 

After a straightforward calculation, one then obtains the formula (15.961 ) stated in the theorem. 
The contour integrals in that formula define a generalized //-function as described in appendix IA1 
The criteria for absolute convergence stated in eq. (IA. 1241) are verified as follows. Let Af be the 
parameter in eq. (IA. 1241) corresponding to the integration variable Wp. We note that \F \ 4>| is the 
number of times that wp appears in the sum Jv,- M&UfbU j) w f » whereas I^V/ 7 ] is the number of 
times that it appears in the sum Y,(ij)e®Y<Ffi(ij) w f- Consequently, e.g. when F G 1e+i 

A F = 2 + 3\F\<S>\-\<i>\F\ =2(l + \F\<f>\) >2, (C.163) 

where in the second step we have used that \F\ = |4>| =V so that |F\4>| = |<J>\F|. In the other case 
when F G % E , we similarly have \F\ = V + 1, |4>| = V and hence \F \ 4>| = |4> \F| + 1, implying 
again 

A F = 1+3|F\4>|-|4>\F| =2(1 + |F\<£|) > 2. (C.164) 
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Therefore, by the criteria stated in appendix |A[ the integrations over the w in eq. (15.961 ) are abso- 
lutely convergent for any |arg(l —Z rs )\ <n< 7linf{A/7 | F G This also justifies aposteriori 
performing the integrations over y in step 5),6) and the changes of integration variables/order of 
integration. 
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